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ABSTRAGT 

We investigate the long-term evolution of X-ray coronae of solar analogs based on high-resolution 
X-ray spectroscopy and photometry with XMM-Newton. Six nearby main-sequence G stars with ages 
between ss 0.1 Gyr and « 1.6 Gyr and rotation periods between « 1 d and 12.4 d have been observed. 
We use the X-ray spectra to derive coronal element abundances of C, N, O, Ne, Mg, Si, S, and Fe 
and the coronal emission measure distribution (EMD). We find that the abundances change from an 
inverse-First Ionization Potential (FIP) distribution in stars with ages around 0.1 Gyr to a solar-type 
FIP distribution in stars at ages of 0.3 Gyr and beyond. This transformation is coincident with a 
steep decline of non-thermal radio emission. The results are in qualitative agreement with a simple 
model in which the stream of electrons in magnetic fields suppresses diffusion of low-FlP ions from the 
chromosphere into the corona. The coronal emission measure distributions show shapes characterized 
by power-laws on each side of the EMD peak. The latter shifts from temperatures of about 10 MK in 
the most rapidly rotating, young stars to temperatures around 4 MK in the oldest target considered 
here. The power-law index on the cooler side of the EMD exceeds expected slopes for static loops, 
with typical values being 1.5-3. We interpret this slope with a model in which the coronal emission 
is due to a superposition of stochastically occurring flares, with an occurrence rate that is distributed 
in radiated energy E as a power-law, dN/dE oc as previously found for solar and stellar flares. 

We obtain the relevant power-law index a from the slope of the high-temperature tail of the EMD. 
Our EMDs indicate a 2.2 — 2.8, in excellent agreement with values previously derived from light 
curves of magnetically active stars. Modulation with time scales reminiscent of flares is found in 
the light curves of all our targets. Several strong flares are also observed. We use our a values to 
simulate light curves and compare them with the observed light curves. We thus derive the range of 
flare energies required to explain the light-curve modulation. More active stars require a larger range 
of flare energies than less active stars within the framework of this simplistic model. In an overall 
scenario, we propose that flaring activity plays a larger role in more active stars. In this model, the 
higher flare rate is responsible both for the higher average coronal temperature and the high coronal 
X-ray luminosity, two parameters that are indeed found to be correlated. 

Subject headings: Stars: activity—stars: coronae—stars: flare—stars: late-type—stars: abundances— 
X-ray: stars 


1. INTRODUCTION 

The solar magnetic field has steadily declined during 
the Sun’s evolution on the main sequence. Studies of stel¬ 
lar clusters and individual field stars with approximately 
known ages have shown that the principal parameter de¬ 
termining magnetic activity on a star is its spin rate 
which, together with convection, controls the operation 
of the internal magnetic dynamo. As a star spins down 
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due to angular momentum loss via a magnetized stellar 
wind, its dynamo action weakens, thus continuously re¬ 
ducing magnetic activity expressed in the outer stellar at¬ 
mosphere. The spin-down history of a solar-like star has 
been studied in detail by using open-c luster samples (e.g., 
lBouvieJII99nHSoderblom et al.j|l 9931) and accompanyii^ 

theoretical calculations (e.g ., iPinsonneault et aLl II989I: 

iMacGregor fc Brenner! 19911) . For a given main-sequence 
stellar mass in the spectral domain of G and K stars, the 
i.diqtMgp^eBiffl^hjS.iahEailhJr well-defined function of the 
star’s age if the latter exceeds a few 100 Myr, regardless 
of the initial Zero-Age Main-Sequence (ZAMS) rotation 
rate. For near-ZAMS stars, in contrast, the rotation rate 
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depends on the pre-main sequence evolution and may 
range between rather modest values (rotation periods P 
of several days) and rates of the so-called ultra-fast ro¬ 
tators (periods of < 1 d), regardless of the precise stellar 
age. This spread in rotation rate is well illustrated by 
stellar cluster samples, such as the Ple iades, in which G 
and K stars are still close to the ZAMS llSoderblom et all 
Il99,11 IStauffer et alJll994jl . 

Magnetic activity expresses itself in the outer stellar at¬ 
mospheres through various phenomena such as dark mag¬ 
netic spots, bright chromospheric plages, chromospheric 
emission lines, and coronal radio and X-ray emissions. 
The coronal emissions display the largest range of vari¬ 
ation in response to the surface magnetic activity level. 
Whereas starspots may cover a few percent of the photo¬ 
sphere of the most active solar analogs, thus producing a 
photometric wave with a modulation depth of a few per¬ 
cent at best, the range of X-ray luminosity Lx in a solar 
analog spans at least three orders of magnitude between 
spun-down inactive examples such as the Sun or (3 Hyi 
(with Lx between a few times 10^® erg s“^ and a few 
times 10^^ erg s“^) and rapid rotators at the “satura¬ 
tion level” (Lx « 10 ~®£hoi » 12 -- 41 x 10®° erg s"^ 
for a solar-mass star; iMagerio et all 1191^ iGiidel et al.1 
n997all . Even the magnetic cycle of the Sun induces an 
X-ray luminosity variat ion over one order of magnitude 
iMicela fc Marinol200.1ll . Such systematic cyclic or semi- 
cyclic variations in turn limit the accuracy with which we 
can attribute “characteristic” X-ray luminosity levels to 
stars of a given age if other solar analogs are subject to 
similar magnetic cycles in X-rays, as seems to be the case 
llDorren fc Guina i1 ll994lborren et alJll99,^ . 

Short-wavelength ultraviolet and X-ray emissions of a 
star like the Sun not only serve as a valuable diagnostic 
to study stellar spin-down and the operation of the in¬ 
ternal dynamo, they are also pivotal for the evolution of 
the outer atmospheres of planets, in particular the pho¬ 
to chemEtr^MuA^ early atmospheres of Ea rth-like plan¬ 
ets llGanuto et all 119^ iRibas et alJl200,^ . Motivated 
by the interest in understanding the Sun’s and our solar 
system’s past, we have been studying the “Sun in Time” 
from the early evolutionary stages on the ZAMS to the 
terminal stage on the main sequence, at ages of 5-10 Gyr. 
This study enc ompasses various wavel ength regimes, in¬ 
cluding r adio llG iidel fc GaidosI 12 00111. optical and ul - 
traviolet llDorren fc Giiinanl 1 1994 iGuinan et ^1200311, 
and X - ray wa velengths llDorren et al.lll995t IGiidel et al.1 
inniiimH). The latter two studies comprise a de¬ 
tailed description of various aspects of the X-ray emis¬ 
sion of a solar analog during its main-sequence evolu¬ 
tion, based on low-resolution X-ray spectroscopic data 
from the ROSAT and ASCA satellites. Among the prin¬ 
cipal findings of that work was a clear trend of decreasing 
coronal electron temperatures as a solar analog ages, fol¬ 
lowing the decrease in overall X-ray luminosity from ages 
of ss 0.1 Gyr all the way to ages of nearly 10 Gyr. The 
authors speculated that the decreasing heating efficiency 
is due to a decreasing level of coronal flaring owing to a 
smaller filling factor in older stars. In this picture, the 
flare rate is responsible both for bringing dense mate¬ 
rial into the corona and heating it to high temperatures, 
analogous to the behavior of individual flares observed 
on the Sun. A large flare rate in younger stars with a 
higher magnetic filling factor could thus produce X-ray 


coronae that are more luminous and at the same time are 
hotter. A larger filling factor is, however, not necessarily 
needed to make the corona more luminous or hotter, as 
we will discuss in this paper. 

The new generation of X-ray observatories, XMM- 
Newton and Chandra, offers entirely new access to stel¬ 
lar coronal physics by providing high-resolution X-ray 
spectroscopy with considerable sensitivity. Apart from 
more detailed studies of the thermal stratification, they 
also allow us to derive the coronal abundances of indi¬ 
vidual elements and to measure electron densities. Both 
parameters are important for our understanding of the 
physical processes in stellar coronae. It is now well 
established that coronae of magnetically active stars 
show various anomalies in their composit ion, in par¬ 
ticular an overall deplet ion of metals (e.g., iDrake et alJ 
IT^ IWhite~^i~^IT^ and a relative underabundance 
of elements with a low (< 10 eV) First Ionization Po¬ 
tential (FIP) com pared to elements wit h a high (> 
10 eV) FIP (e.^. IBrinkman et aP 120011: IGiidel et alJ 
l2001llDrake et alJ2001ll . In the inactive Sun, in contrast, 
element abundances are arranged according to the so- 
called FlP-effect, with low-FIP elements being overabun¬ 
dant in the corona relative to their pho tospheric values 
and relative to the higher-FIP elements llMeveilll985albl : 

m Simil ar trends have been noted in in¬ 

active stellar coronae llTarniug et al 1119961: IDrake et all 
IT^lGiidel et alJ2002ll . Because the elemental composi¬ 
tion and the element fractionation ultimately derive from 
the photospheric gas and the physical processes heating 
it and transporting it to coronal heights, they may be im¬ 
portant tracers for the coronal heating mechanism. For 
example, it has been suggested that the Ne enrichment 
seen in active stars may be re lated to increased levels of 
flaring llBrinkman et al.ll2001ll . 

This paper presents a study of six solar analogs at 
young and intermediate ages observed with the instru¬ 
ments onboard XMM-Newton. The stars range from 
near-ZAMS ages (0.1 Gyr) to ages between approxi¬ 
mately 1-2 Gyr when the rotation period has increased 
to about 12 days. Together with the Sun, they cover al¬ 
most the entire path of coronal main-sequence evolution. 
The present paper is in many ways a follow-up and con - 
tinuation of the work presented bv IGiidel et al.l (^223) ■ 

Our work emphasizes the similarity of stellar mass, 
surface gravities and internal structure, i.e., we confine 
this study to the main sequence and treat rotation (or 
equivalently, age) as ou r principal free par ameter. In 
a complementary study, iScelsi et"^ ll2005ll investigate 
three G-type stars, all at very high activity levels but 
at largely differing evolutionary stages and with differ¬ 
ent internal structure and surface gravities, including the 
young main-sequence star EK Dra, a pre-main sequence 
(weak-lined T Tau) star, and the evolved Hertzsprung- 
gap giant 31 Com. They find very similar emission mea¬ 
sure distributions in the latter two active stars, regardless 
of the differences in their fundamental parameters. 

The outline of our presentation is as follows. In Sect. [3 
and 13 we discuss our targets, the observations, and the 
principal procedures of the data reduction, respectively. 
Special attention has been paid to the analysis and inter¬ 
pretation of the spectroscopic data in order to recognize 
the principal potential and the limitations of the spec¬ 
tral inversion, i.e., the derivation of emission measure 





























































Coronal Evolution of the Sun in Time 


3 


distributions and the element abundances. In Sect. 0 we 
describe the spectra of the six targets. We discuss our 
methods in Sect.El Sect. [HI presents our results, while we 
discuss various features and models in Sect. [7| Finally, 
Sect.lHl contains our conclusions. 

2. TARGETS 
2.1. General Properties 

XMM-Newton data of six young and intermediate- 
age solar analog stars have been analyzed. The stars 
are all of early-to-mid spectral type G on the main se¬ 
quence. Their ages range from approximately 0.1 Gyr 
for 47 Cas B and EK Dra to « 1.6 Gyr for [3 Gom. 
These ages have been determined using various prox¬ 
ies such as rotation periods (for the older targets), or 
memberships in moving groups of known ages (for the 
younger targets). We have selected these targets be- 
cause they have been st udied in much detail before (e.g., 
iDorren fc Guinai]ll994D . have well-measured fundamen¬ 
tal parameters, and are well-behaved representatives of 
their age cl ass. A detailed s umma ry an d discussion 
is give n in IDorren fc GuinaiJ lll994H and iGiidel et al.1 
The dist ances quoted are derive d from the Hip- 
parcos parallaxes ijPerrvma.n et a,]]l1997jl . and the rota¬ 
tion periods have been measured photometrically. The 
general properties of the stars are listed in Table Q] 
where they are also compared with solar values. The 
X-ray luminosities Lx given there refer to measure- 
ments from ROSAT photom etry in the 0.1-2.4 keV band 
iGiidel et al.ll997^ll998all . and to the spectral modeling 
discussed in the present paper, also for the 0.1-2.4 keV 
band, and additionally for the 0.1-10 keV band. For 
the X-ray luminosity of the S un, we use a represe ntative 
value of log Lx = 27.3 as in IGiidel et al.l l|1997ali . This 
val ue is in agreement wi th the Lx of a Gen determined 
bv lR,aassen et al.l (j2003), Lx = 1-6 x 10^^ erg s“^. A 
steady decline of Lx with increasing age and increasing 
rotation period P is evident. 

The most active target, 47 Gas B, is the fainter com¬ 
ponent in the 47 Gas binary system. It has not been 
characterized optically as it h as been individually de¬ 
tected only by radio methods (IGiidel et alJll998bfi and 
indirectly from Hipparcos measurements as a compan¬ 
ion to the optically bright FO V star 47 Gas = HR 581 
= HD 12230. The radio position is clearly offset from 
the position of the FO star. It is, furthermore, very un¬ 
likely that early F st ars display luminous and spectrally 
hard X-ray emission (iPanzera et a,1 .1(^19911 . On the other 
hand, all characteristics of the X-ray detection fit well to 
an early G-type star with an age similar to that of the 
Pleiades. In particular, the X-ray luminosity corresponds 
to the saturation level of an early G star (Lx « 10“^Lboi, 
lyiThllT^ . thus excluding a corona of a later-type star 
as its origin. 

A saturated corona requires, for a G star, a rotation 
period smaller than that of EK Dra. A periodic signal 
with a period of w 1 day was reported from the ROSAT 
All-Sky Surv ey observations, an d was attributed to stel¬ 
lar rotation (IGiidel et abllT99^ . This period also coin¬ 
cides with the rotation p eriod of the fastest earl y-G type 
rotators in the Pleiades ()Soderblom et al.lll993li. We re¬ 
fer th e rea der to the de t ailed d iscussion in lGiidel~ iiO 
(I1995^ and IGiidel et alJ (Il998bli . We tried to quantify 
the maximum contribution of the FO star to the X-ray 


luniinosit |^ of the 47 Gas system. ^ ROSAT studi^ of t^ 

found only few stars with spectral type A7-F3 that have 
log Lx > 29.4 (i.e., more than 10% of the Lx of 47 Gas), 
and all of them have known late-type companions, which 
are likely to be the sour ces of the high X-ray luminosity 
(iMermilliod et al1ll992H . In addition, the typical X-ray 
spectrum of F-type s tars in the Pleiade s is softer than 
that of G-type stars (jGagne et al.lll99,^ . so we expect 
the F primary to provide even less flux from hot plasma. 
In summary, we thus expect the FO star to contribute less 
than 10% to the total detected X-ray flux of the 47 Gas 
system. We therefore add this target to our list of G-type 
stars, being the only likely solar analog that is accessible 
to high-resolution X-ray spectroscopy and that is more 
active than EK Dra, in fact reaching the saturation level. 

2.2. Photospheric Composition 

When measuring coronal abundances, A, of elements, a 
principal problem is to what standard set of abundances 
the results should refer to. Often, stellar coronal abun¬ 
dances are cited with respect to the solar photospheric 
composition. But because the stellar coronal plasma ul¬ 
timately originates in the respective stellar photosphere, 
coronal abundances may reflect the composition of the 
latter, and abundance “anomalies” may simply be due 
to a non-solar composition of the stellar photospheric 
material. Fortunately, the photospheric abundances of 
most of our objects have been measured. The reported 
photospheric abundances of our sample of stars (except 
47 Gas) are given in Table|21 The entries refer to individ¬ 
ual measurements or to catalogs compiled from previous 
studies. 

As one sees from this literature survey, almost all Fe 
abundances reported for our targets are compatible with 
solar photospheric values. Moreover, several measure¬ 
ments for other elements exist, including high-FIP ele¬ 
ments such as G, N, and O, and again most of the re¬ 
ports are compatible with solar abundances. We find 
a trend for super-solar abundances in [3 Gom, by per¬ 
haps 10-20%. Nevertheless, the photospheric-abundance 
tabulations for our targets are incomplete and the given 
measurements scatter, making it impossible to reliably 
express our coronal abundances relative to the respec¬ 
tive photospheric values. However, the above summary 
makes it clear that our targets (except for 47 Gas) must 
show a near-solar composition, and we therefore adopt 
the solar photospheric abundances as the reference com¬ 
position throughout the paper. As for 47 Gas, to our 
knowledge no photospheric abundances of the brighter 
F star in the system have been reported. However, this 
syste m is young and is a member of the Local Associa¬ 
tion (IGiidel et al.lll998bh . for which we can reasonably 
assume near-solar abundances, similar to EK Dra, which 
is also a member of the Local Association. 

3. OBSERVATIONS AND DATA REDUCTION 
3.1. Observations 

Our target stars were obser ved with the Ref l ection 
Grating Spectrometers (RGS, Iden Herder et al iunni 
and the European Photon Imaging Gameras (EPIG, 
iStriider e t alJl200lHTurn er et alJl2001ll onboard XMM- 
Newton (iJansen et al .1 120011) . The RGS are suited for 
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high-resolution spectroscopy in the wavelength range be¬ 
tween 6-38 A, with a resolution of A A « 60 — 76 mA, 
hence a resolving power of A/AA ~ 100—500. The EPICs 
observe between « 0.15-15 keV, providing a moderate 
energy resolution of approximately E/AE = 20 — 50. 

We provide a detailed log of the observations in Ta¬ 
ble 13 

3.2. Data Reduction 

The data were reduced using the Science Analysis Sys¬ 
tem (SAS) version 5.4.1. We applied the standard pro¬ 
cessing performed by the RGS metatask rgsproc and 
the EPIC MOS task emchain. The PN data (used only 
for the light curves) were reduced using the PN task 
epchain. 

For the RGS data reduction, we extracted the 
first-order spectra from a spatial cut including 95% 
of the cross-dispersion Point Spread Function (PSF) 
(xpsfincl=95 in rgsproc) and an energy cut includ¬ 
ing 95% of the pulse-height distribution (pdistincl=95). 
The background spectra were extracted above and below 
the source spectra, by excluding 97% of the source cross¬ 
dispersion PSF. 

For all stars except UMa, O^i, and (3 Com the 
MOS data were taken in the small-window mode. In this 
mode, only the central 100 x 100 pixels of the middle CCD 
are operational, whereas the outer CCDs work in full 
mode. The small-window mode allows for shorter inte¬ 
gration times for the source and thus avoids the brighter 
sources from becoming piled up. The disadvantage of this 
mode is, however, that the central window is too small 
to reliably extract a background region. We therefore 
selected a source-free region in one of the outer CCDs. 
The target source itself was extracted from a circle as 
large as possible in the small window (with a radius of 
approximately 50"). 

For TT^ UMa, MOSl was in full window mode whereas 
MOS2 was in small-window mode. The exposure time 
was 5000 s longer in MOSl. Pile-up was studied for both 
MOS cameras with the SAS epatplot task. This task 
makes use of the relative ratio of single- and double-pixel 
events, that deviate in case of pile-up. Since no evidence 
for pile-up was found for this star, we decided to use 
MOSl. 

For x^ Ori, MOSl data were taken in timing mode 
and MOS2 data in full-window mode. We did not use 
data taken in timing mode because reliable background 
subtraction cannot be performed. We found that pile-up 
was not negligible in the MOS2 data and therefore used 
an annular extraction region for the source, with inner 
and outer radius of 80 and 1200 detector pixels (4 and 
60 arcsec), respectively, in order to remove the piled-up 
central part of the PSF. 

For P Com both MOSl and MOS2 observed in full 
window mode. We found that pile-up did not affect this 
source. 

In those cases where the background area was se¬ 
lected in one of the outer CCDs, the background count 
rates could have been underestimated with respect to the 
source extraction region due to vignetting. We studied 
background behavior effects using the eboxdetect SAS 
task. This task detects all sources in the field of view; 
these were subsequently cut out from the image. We then 
quantified the remaining background count rate per unit 


area as a function of the distance from the center of the 
field of view, by analyzing the average rates in annuli 
at different radii. We found a decrease of at most 10% 
from the image center to the central portions of the outer 
CCDs, which is negligible when compared to calibration 
uncertainties, particularly in the light of the high count 
rates of our sources. The only possible but small influ¬ 
ence could be on the data at the highest energies where 
the source flux drops below the background flux. We 
therefore decided not to consider spectral data at ener¬ 
gies where the background count rate exceeded the source 
count rate. 

We used the rgsrmfgen task for RGS and the rmfgen 
and arfgen tasks for MOS to generate response and an¬ 
cillary files appropriate for the specified extraction re¬ 
gions. To avoid bins with very few counts, we grouped 
the spectral channels: the MOS data were grouped to a 
minimum of 25 counts per bin, whereas the RGS data 
were grouped to a minimum of 10 counts per bin. The 
very weak continuum level in the RGS spectra of P Com 
made this grouping scheme not well suited; instead, we 
regularly rebinned the entire RGS spectra by a factor of 
five. 

Finally, we studied the agreement between the ob¬ 
served wavelengths of bright emission lines with their 
tabulated rest wavelengths (stellar radial velocities be¬ 
ing negligible); in cases where small systematic offsets 
due to residual calibration uncertainties of order 10 mA 
were found, we slightly updated the assumed boresight 
coordinates accordingly and repeated the data reduction. 

4. DESCRIPTION OF THE SPECTRA 

Before quantitatively evaluating the observed spec¬ 
tra, we discuss some general features. Fig. ^ re¬ 
veals four significantly different types of spectra: 
i) the most active star, 47 Cas B, shows typi¬ 
cal features of a very hot plasma, namely a well- 
developed bremsstrahlung continuum, lines of Mg XI 
and XII, and comparatively high flux ratios of 
Oviii A18.97/OVII A21.60, of Nex A12.13/Neix A13.45, 
and of Fexviii A14.20/Fexvii A15.01. ii) The spec¬ 
trum of EK Dra is significantly cooler, which is in par¬ 
ticular evident from the smaller Oviii/Ovii flux ratio, 
a more modest continuum level compared to the line 
strengths, and a marked dominance of the Eexvii lines, 
iii) The third group consists of the intermediately ac¬ 
tive, intermediately old stars tP UMa, x} Ori, and 
Get which all show very similar spectra in which the 
Fexvii lines are much stronger than those of Eexviii, 
and the Oviii/Ovii flux ratio is further reduced. Note 
that the maximum formation temperature of EeXVii is 
only « 5 MK. The continuum has become very weak 
in these spectra, iv) Finally, the spectrum of the least 
active target, P Com, is largely dominated by lines of 
FeXVii, and the O Viii/0 vii flux ratio approaches unity. 

Turning to the EPIC MOS spectra in Fig. [3 further 
indicators support this picture. In this figure, we have 
renormalized the MOS spectra so as to represent the stars 
at a common distance of 33.56 pc, identical to the dis¬ 
tance of 47 Cas B. A marked decrease of the overall emis¬ 
sion level is seen as the stellar age proceeds from the sam¬ 
ple of Pleiades-age stars (47 Cas B, EK Dra) to the older 
sample. The most active stars show shallower continuum 
spectra between « 2 and 10 keV than the less active tar- 
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gets, indicative of the higher overall plasma temperatures 
of the former. The dominant Fexvii lines (at 0.826 keV 
and 0.727 keV) are well developed in the spectra of the 
less active stars. 


5. DATA ANALYSIS 

Because the novel aspect of our X-ray data is the high 
spectral resolution of the RGS, allowing us to access in¬ 
dividual emission lines from different elements, we used 
the spectral data from both RGS instruments but used 
EPIC data only in so far as they contribute additional in¬ 
formation not accessible by RGS, such as spectral data at 
wavelengths shorter than 6 A. To keep maximum weight 
in our data analysis on the RGS data, we restricted the 
more sensitive EPIC spectral data to one of the three 
cameras. We chose MOSl or MOS2, because the MOS 
spectral resolution is superior to that of the pn camera, 
and the S/N at higher energies balances well with the 
S/N provided by RGS at lower energies. Also, the cross¬ 
calibration between MOS and RGS is better understood 
(they use the same mirrors). Only for the light curves 
did we make use of the data of other EPIG cameras as 
well. 

For the analysis of each target, we used the exposure 
time covered simultaneously by all three detectors. Each 
light curve except that of (3 Gom contained one well- 
developed flare. In order to avoid systematic bias by 
any of these flares (i.e., increasing the average Lx and 
possibly also increasing the characteristic coronal tem¬ 
peratures), we eliminated the flare intervals from consid¬ 
eration. We will, however, briefly present and discuss the 
individual light curves separately in Sect. lO The expo¬ 
sure times remaining for our spectral analysis are listed 
in Table 0 (last column). 

5.1. Spectral Inversion 

An observed coronal X-ray spectrum is the superposi¬ 
tion of the spectra emitted by various coronal features 
with different temperatures, volumes, densities, and pos¬ 
sibly different composition in terms of chemical elements. 
The inversion of an observed spectrum to obtain the un¬ 
derlying physical parameters is therefore a highly degen¬ 
erate problem, with numerous solutions describing essen¬ 
tially the same spectrum. 

For our analysis we consider a coronal model with 
the following, observationally tested assumptions of the 
physical parameters. First, the plasma is assumed to 
be in collisional ionization equilibrium, a model that ap¬ 
pears to be sufficiently good as long as no very rapid 
change in th e heating rate of the plasma is taking place 
l|Mewelll999jl . We also assume that the coronal plasma 
is effectively optically thin and that lines with high oscil¬ 
lator strengths are not subject to resonance scattering, 
an assumption that has been shown to be justified for 
the stars in our sample ijNess et 3.1.1120(1,'^ . Furthermore, 
the density-dependence of the populations of metastable 
levels is neglected, i.e., the spectrum is computed in the 
low-density limit. This approximation is supported by 
the flux ratios in the He-like triplet of O vii in our ob¬ 
servations, and it appears t o be a reasonable assumption 
for m ost coronal plasmas l|Ne.ss et a,1.ll2flf)4l iTesta. et iiD 

l2finl . 

At this point, then, the observed spectrum is essen¬ 
tially a function of the distribution of optically thin coro¬ 


nal features in volume, temperature, and elemental com¬ 
position; conventionally, the thermal distribution is de¬ 
scribed by the differential emission measure distribution 
(DEM), 


Q{T) 


UgUndV 

dlnT 


( 1 ) 


where Ue and uh are the electron and hydrogen number 
densities, respectively, and dV is a differential volume 
element at temperature T. The DEM determines the 
line flux Fj of any given emission line j through 


A<p,{T)QiT)d\nT. (2) 


Here, ipj is the line power per unit emission measure 
(“emissivity” henceforth), and A is the abundance of the 
element emitting the relevant line, with respect to some 
standard tabulation (such as “solar photospheric abun¬ 
dances” ) used for the computation of (pj. 

Non-solar abunda nces in stellar coronae are now 
well established ('e.g..lBrinkman et al.ll200lHDrake et alJ 

l200HlATidard et alJ2003all . and we recall that the corona 

of the Sun itself shows considerable deviations from the 
photospheric mixture. However, it is also known that 
the solar coronal abundances vary greatly from feature 
to feature, some showing a marked FIP bias, with oth¬ 
ers showing photo spheric composition llFeldmanI 119921: 
iFa.ming et ablfl 99511 . Our spectra are - like most existing 
stellar X-ray spectra - insufficient to characterize abun¬ 
dances at various separate temperatures; therefore, we 
will assume T-independent abundances to at least recog¬ 
nize overall trends in our stellar sample. 

The spectral inversion problem is mathematically ill- 
posed. Statistical scatter due to photon statistics in line 
fluxes, even if amounting to only a few percent, may 
introduce considerable sca tter in the reconstruc ted emis¬ 
sion measure distribution llGraig fc Brownll976ll . We ex¬ 
pect a strong amplification of such effects from the un¬ 
certainties in the atomic physics parameters (which may 
amount up to a few tens of percent for the emissivities) 
and from calibration uncertainties (up to several percent 
in certain spectral regions). 

Therefore, the spectral inversion can essentially be 
performed in a meaningful way only by constraining 
the problem suitably, for example by subjecting the re¬ 
construction to smoothness conditions for the result¬ 
ing, discrete Emission Measure Distribution (EMD), or 
by iterating the problem using a pre-defined conver¬ 
gence algorithm that reconstru cts a preferred, physic ally 
meaningful solution (see, e.g.. iKaastra et af]ll996al and 
iGiidel et alJll997bl for a discussion of various algorithms 
with different bias). 

We have chosen here to perform the inversion using two 
widely different algorithms that we briefly summarize as 
follows: 


• Method 1 (Ml). We fitted the spectrum using syn¬ 
thetic template spectra calculated for a set of phys¬ 
ical parameters; the parameters were varied until 
the fit was optimized. This method has conven¬ 
tionally been used for low-resolution spectra but 
also for high-resolution grating spe ctra, in partic¬ 
ular a s implemented in t he SPEX llKaa stra, et all 
Il996bll and the XSPEC llArnaudI 1199611 software 
packages. We, however, modified the conventional 
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approach by including almost exclusively segments 
of the spectrum that are dominated by bright lines 
for which comparatively robust atomic physical pa¬ 
rameters should be available, together with some 
nearly line free regions. 

• Method 2 (M2). Here, we worked with a list of 
discrete line fluxes that we extracted from the ob¬ 
served spectrum. If the formation temperatures of 
these lines cover the range of relevant coronal elec¬ 
tron temperatures, then the EMD can be recon¬ 
structed by successive approximation, essentially 
inverting a system of equations like Eq. 0 . 

We emphasize the complementarity of our two ap¬ 
proaches: method 1 uses all tabulated emission lines and 
their blends within the selected spectral segments, while 
we confine method 2 to a minimum number of lines, 
namely the brightest lines available, required to derive 
meaningful EMDs and the most essential abundances. 

Both methods are subject to separate biases. Some of 
the many lines considered in method 1 may be poorly fit¬ 
ted, owing to discrepancy in the tabulated emissivities. 
If the latter scatter around their true values, the result¬ 
ing fit may show reduced systematic errors compared to 
a fit based on one single line per ionization stage. The in¬ 
dividual lines in method 2 may be better described, but 
their small number makes the inversion process rather 
sensitive to any systematic uncertainty in any of the line 
emissivities, and some a priori estimate for blend contri¬ 
bution must be considered. Lines not tabulated in the 
line-emissivity lists affect both methods, either by alter¬ 
ing the measured line flux in method 2 through unrecog¬ 
nized blending, or by adding a continuum-like base level 
if many weak lines contribute in either of the analysis 
methods. 

A comparison of method 2 with a method using poly¬ 
nomial EMDs to calculat e synthetic spectra h as previ¬ 
ously been presented by lAudard et al.l l|2004H for the 
extremely active and hot corona of the FK Corn-type 
star YY Men. These authors found that both the re¬ 
constructed EMDs and the derived abundances closely 
agreed for the two methods. Here, we study and com¬ 
pare the results for stars across a large range of activity 
and coronal temperatures. 

To further assess potential sources of systematic error, 
we applied two different databases to each of the two 
methods. The spectral models are based on the assump¬ 
tion of coronal ionization equilibrium (CIE). We used the 
atomic p arameters from the MEKAL database in S PEX 
vers.2.0 llMewe et al.l 119851: iKaastra et al.l Il996bll and 
from the Astrophysical Pla sma Emission Code vers. 1.3. 1 
fAPEC. iSmith et abl 1200111 in XSPEC (|Arna,iid| 11 QOOfl . 
The MEKAL emis sivities are computed using the ion- 
izat ion balance of [Arnaud fc R avmondl lll992H for Fe 
and lArnaud fc R.othenflusI l|1 98.511 for the other elements. 
The APED emissivities (emissivities used in the APEC 
code) are computed using the ionization balance of 
iMazzotta et al.1 i|1998ll . 

5.2. Method 1: Synthetic Spectra 

Here, we used an approach essentially identical to 
the one described in lAudard et al.l ll2003ar . apart from 
slightly different wavelength ranges (see below). Because 


numerous emission lines are poorly described by the spec¬ 
tral line lists in use, we confined our analysis to a re¬ 
stricted number of spectral regions that contain bright 
lines for which the atomic physics is believed to be rela¬ 
tively well known, and some nearly line-free regions of the 
spectrum where the continuum dominates. The wave¬ 
length ranges of the selected regions are summarized in 
Table m As we decided to put more weight on the data 
from the RGS with its higher spectral resolution, we used 
EPIC MOS data only in the wavelength range shortward 
of 9.35 A, where the RGS effective area is small and its 
calibration is less accurate. We thus cut the RGS data 
shortward of 8.3 A, so that the data from the RGS and 
the EPIG instruments overlapped around the Mg XI and 
Mg XII lines. As mentioned in Section 4, we discarded 
the high-energy part of the EPIG spectrum where the 
background flux exceeded the source flux (see Table El). 

The physical model was defined as a combination of 
10 optically thin, thermal CIE models, and a photoelec¬ 
tric absorption component. The photoelectric absorption 
was frozen by defining a fixed interstellar hydrogen col¬ 
umn density. The latter ranged between negligible val¬ 
ues (< 10^® cm“^) for the closer stars and 7 x 10^® cm“^ 
for 47 Cas B. These values a re consistent with hy dro¬ 
gen column d e nsitie s given bv lGiidel et al.1 lll997ljl and 
lAudard et HI lll999ll and with va l ues of the interstellar 
hydrogen density given byHirigiilEnil). However, even 
the largest values used here were too small to significantly 
alter our model spectra even at the longest wavelengths 
in the RGS. 

Each of the ten GIE temperatures was confined to 
within one of ten adjacent temperature intervals of 
equal width (0.22 dex) in logT; to optimize the multi¬ 
temperature fit, given the relatively broad bins, we treat 
the temperature values as free parameters within the 
bounds of the respective bin intervals, as well as the as¬ 
sociated emission measures. The temperature intervals 
covered the range from logT=6 to logT=8.2 (T is given 
in K). The abundances of G, N, O, Ne, Mg, Ar, Si, and S 
were also used as free parameters. Only the abundances 
found with the 10-temperature fit were used in a second 
step for the DEM reconstruction. The Ghebychev poly¬ 
nomial DEM code in SPEX, with polynomials of order 
6 and 8, was used to describe the DEM that best fitted 
the given spectrum. An example of a best fit is shown in 
Fig.El which also illustrates the spectral ranges used for 
this procedure. 

5.3. Method 2: Inversion of Line-Flux Lists 
5.3.1. EMD Reconstruction 

We reconstructed the discrete EMD starting from 
fluxes of individual lines. In order to obtain an EMD 
independently of the element abundances, we selected 
a few well-defined, bright Fe lines from the ionization 
stages of Fexvii to Fexxv (as far as measurable, see 
Table EJ . The Fe xvii-xxiv lines were extracted from 
the RGS spectrum, whereas the unresolved Fexxv line 
triplet was extracted as a single blend system from the 
MOS spectrum. The emissivities of the Fexvii-xxiv 
lines cover the temperature range between 3 MK and 
20 MK rather well, and for the more active stars (47 Gas 
B and EK Dra), Fexxv constrains the hottest tempera¬ 
tures, up to about 100 MK (Figure^. In order to obtain 
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information on the cooler part of the EMD in which no 
Fe lines are detected, we used the flux ratio between the 
OVIII A18.97 and the OVII A21.60 resonance lines, which 
itself is independent of abundances. The lines used for 
each star are listed in TableEland TableEl together with 
their measured luminosities. 

To avoid cross-talk with the abundance determination, 
we selected only lines of Fe that are not strongly blended 
with lines of other elements (with the exception of Fe xix, 
see below; the Fexvii lines at 17A were not used be¬ 
cause they are partly obscured by one of the RGS CCD 
gaps). Nevertheless, our line features usually still contain 
a number of weaker blends from various Fe ions. We thus 
considered all Fe blending lines close to the principal line 
as being part of our Fe feature and we therefore computed 
new T-dependent emissivity curves for each such Fe hlend 
system. The line feature, composed of all Fe lines in the 
chosen wavelength range, was thus subsequently treated 
like a single emission line. In the RGS spectrum, all tab¬ 
ulated blending Fe lines within ±0.06 A (one FWHM on 
each side) of the wavelength of the principal line were 
considered. In the MOS spectrum, the range was chosen 
to be ±0.15 A. 

In a first step, we extracted the fluxes of the relevant 
Fe lines and of the Oviii A18.97 and Ovii A21.60 res¬ 
onance lines from the spectrum. This can be problem¬ 
atic, because the lines have broad wings due to the RGS 
PSF, and the determination of the underlying contin¬ 
uum is difficult. Moreover, in some cases, blends with 
lines from different elements may still be present, in par¬ 
ticular in the case of the brighter lines of Fexix that 
are blended with Neix lines. While the effect of these 
blends cannot be assessed a priori without knowledge 
of the thermal structure and the abundances, we can 
obtain useful approximations as follows. Because most 
ions are represented by several different emission lines, 
and the continuum can be interpolated from values ob¬ 
tained in nearly line-free regions of the spectrum, a good 
approximation of the contaminating blends and of the 
baseline continuum level in fact comes from the fit we 
derived in method 1. We emphasize that we have not 
used that model to derive parameters, but only to ob¬ 
tain an approximate description of the baseline spectral 
flux distribution on which the line fluxes of interest are 
superimposed. The latter being the fluxes of Fe lines, we 
set the Fe abundances in the best-fit model from method 
1 to zero, retaining only the emission lines and the con¬ 
tinuum from all other elements. We note that neglecting 
Fe does slightly change the continuum as well, but the 
influence is on the order of a few percent in the tempera¬ 
ture and wavelength range of interest if solar abundances 
are assumed. Anticipating sub-solar abundances in ac¬ 
tive stars with a strong continuum as reported in the 
previous literature, the contribution of Fe to the con¬ 
tinuum will be even smaller. The errors thus introduced 
into the line-flux extraction are negligible. An equivalent 
procedure was then also applied to extract the O viii and 
O VII line fluxes, although in these cases, blending was 
not significant. 

Before extracting individual lines, we need to make 
sure that the baseline continuum level has been correctly 
fitted individually in each spectrum. We checked the 
relevant fits in a nearly line-free region of the RGS spec¬ 


trum, namely immediately shortward and longward of 
the Oviii A18.97 resonance line by multiplying the EM 
by a suitable factor. We found that an optimization of 
the continuum level in these wavelength regions required 
a change in flux of only a few percent. It is not clear 
whether this is due to slight cross-calibration problems 
or to some bias in the fit. Such bias may be introduced 
by the fact that the MOS spectral portion, determin¬ 
ing mostly the hotter part of the DEM, has a relatively 
high S/N ratio which may lead to some over- or under¬ 
estimation of the hot EM. Because the latter produces 
continuum at all wavelengths, a slight offset could also 
affect the soft part of the RGS spectrum. Given the small 
magnitude of the effect, we are unable to ascribe such an 
offset either to (possibly wavelength-dependent) cross¬ 
calibration problems or to fit bias as described above. 
Our final results indicate that any such effect is of minor 
consequence and does not need to be considered for the 
overall fit in method 1. For the extraction of line fluxes 
in method 2, however, we need to optimize the contin¬ 
uum level so that the line flux can be properly defined 
as an excess above this baseline level. 

In a similar manner, we extracted the Fe xxv A1.85 line 
by first optimizing the nearby continuum level starting 
from the overall best-fit provided by method 1. Again, 
the continuum renormalization required for this proce¬ 
dure was at most a few percent in the spectral region of 
interest. 

To actually extract the line fluxes, we added 5 function 
lines to the model, located at the theoretical wavelengths 
of the principal lines in the blends under consideration, 
and convolved them with the instrument response. The 
line fluxes were then obtained by fitting the amplitudes 
of the S functions to the observed lines and the adja¬ 
cent narrow wavelength ranges. If the wings of two Fe 
lines overlapped significantly, we fitted two 6 functions 
simultaneously. Note that the fluxes thus derived may 
be different for APEC and MEKAL because of slightly 
different atomic physics, fit parameters, and blend con¬ 
tributions. 

The obtained Fe, O viii and O Vii line fluxes were con¬ 
verted to luminosities as presented in Table 0 and Elfor 
MEKAL and APEC emissivities, respectively, based on 
Hipparcos distances as quoted in Table Q They were 
then used to reconstruct the EMD, as follows: 

A first, smooth estimate of the EMD was derived from 
the emissivities at the maximum line formation temper¬ 
ature Tm of each Fe line. T he emissivities are based o n 
solar abundances as given bv I Anders fc Grevessd l|1989ll . 
Our EMD was defined on a grid with a grid point sepa¬ 
ration of Alogr=0.1 dex in the range between logT = 
5.5 - 8 for SPEX (T in K) (APEC: logT = 5.7 - 8) for 
the more active stars 47 Gas B and EK Dra, and in the 
range between logT = 5.5—7.1 (APEC: logT = 5.7—7.1) 
for the less active stars Ori, Get and (3 Com. For 
TT^ UMa, the EMD was defined in the range between 
logT = 5.5 - 7.3 (APEC: logT = 5.7 - 7.3). (The dif¬ 
ferent low-temperature limits reflect the different avail¬ 
ability of emissivities in the different codes, but we note 
that the range of logT = 5.5 — 5.7 is irrelevant for any 
line we use for the EMD reconstruction, and it is also 
not important for the emissivities of any lines of other 
elements analyzed here.) As a starting condition, the 
EMD was extrapolated to temperatures cooler than Tm 
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of Fexvii (logT = 6.7) by using a slope of 2, and to 
temperatures higher than Tm of the hottest Fe ion in use 
(Fexxv, logT = 7.9 for 47 Cas and EK Dra, Fexix, 
logT = 6.9 for P Com, and Fexxi, logT = 7.0 for the 
other stars) by using a slope of -2. These extrapolations 
were performed to the limits of the respective tempera¬ 
ture ranges as defined above. These starting conditions 
were suggested by the slopes of the initial EMDs defined 
by the EM values at the different T^. Also, we subse¬ 
quently found that the low-T slope converges to values 
around 2 even if the starting slope was largely different 
(see Sect. rot . Eor each line, the flux Fq predicted from 
the initial EMD was calculated according to the following 
equation: 

Fc{x) ='^EM{i) ■ ip{x,i) (3) 


where the sum is over the temperature bins i; EM(i) is 
the emission measure in the ith temperature bin, and 
(p{x, i) is the emissivity of the given line x at this grid 
point. The calculated fluxes were then compared with 
the measured fluxes, F{x). The EM in each bin was 
iteratively corrected using the algorithm described by 
iWithbrofi iTmTl : 


EM”+i(f) = EM”(ii)- 


E, 


Y.xV{x,i) -b<p(OVII,i) 


( 4 ) 

where i is the index of the temperature bins and n is 
the iteration number. The last term in the numerator 
and denominator was added to also obtain convergence 
of the line-flux ratio between O viii and O vii, where we 
have used the measured (i?) and the predicted {Rc) flux 
ratio. The EMD was then iterated until we reached a 
pre-set convergence criterion. 

In the presence of considerable systematic uncertain¬ 
ties in the modeling, a reasonable convergence criterion 
should be set, although its statistical meaning may be 
marginal. A straightforward goal is to achieve conver¬ 
gence in such a way that on average the squared devia¬ 
tion between model and fit for a line flux is equal to the 
variance of the same line flux. This leads to a reduced- 
X^-like expression. 


{F[x]-Fc[xm {R-Rcr \ 

(^N)^ J ) 

( 5 ) 

where the last term relates to the Oviii/Ovii flux ra¬ 
tio. Here, is the number of lines used, and a{x) is 
the error assigned to the measured line flux, which com¬ 
prises the error from the finite photon statistics and an 
assumed systematic uncertainty from the atomic physics 
(see Sect. 15.3l^ . Similarly, an is the derived error in the 
Oviii/Ovii flux ratio. We iterated until this expression 
reached a value of unity, or if this did not occur, until it 
no longer significantly decreased. We would like to men¬ 
tion that the fit parameters, i.e., the EMs per bin, are 
not independent of each other owing to the broad emis¬ 
sivity curves of each line; given the large systematic un¬ 
certainties in the line fluxes, we feel that a more detailed 
convergence criterion is not warranted. Our main goal is 
to stop the iteration at a reasonable level to avoid over¬ 
interpretation of superhcial features in the EMD that 




1 


Nx + 1 


E 


may arise from iterating too deeply - see our discussion 
in Sect.O 

At this point, then, we have found an EMD from fluxes 
of the single element Ee, under the assumption of solar 
metallicity as used for the tabulated emissivities, and 
from a ratio of O line-fluxes. Therefore, there still re¬ 
mains a normalization factor for the EMD and the ab¬ 
solute level of the Fe abundance to be determined (see 
Sect.ES3- We will use the observed continuum level to 
fix the EMD normalization and, at the same time, the 
absolute Ee abundance, as explained in the following. 

5.3.2. Abundances 

To determine the abundances of the elements, we ex¬ 
tracted all lines of interest (see Tables 0 and 0 in a sim¬ 
ilar manner as we extracted the Fe and O lines before. 
Here, however we used the EMD constructed from our Ee 
and O line fluxes in order to describe the continuum, in¬ 
stead of the 10-T model. This EMD is known only up to 
a normalization constant depending on the absolute Fe 
abundance, which we derived as follows: we constructed 
a set of spectra from the calculated EMD with different 
Fe abundances such that the product of the Fe abun¬ 
dance and the EM(T) is constant for any T, still using 
the approximate abundances of the other elements from 
method 1 to estimate the contributions from blends and 
to obtain a more accurate description of the continuum. 
From the spectrum that best fitted the nearly line-free 
regions long- and shortward of the Oviii Lya line at 
18.97 A, we obtained the absolute Fe abundance. We 
then extracted the line fluxes of all interesting elements 
using (5-line models, as described before. Note that the 
line fluxes of Mg, Si, S, and, if available, Fexxv, were ex¬ 
tracted from the MOS spectra while all other lines were 
extracted from the RGS spectra. As done for the RGS 
spectra before, we adjusted the model continuum to the 
observation in the line-free regions of the MOS instru¬ 
ment at high energies, before the line-flux measurement. 

The predicted line fluxes of these elements (see TableEl 
and E|) were then calculated from the EMD, using their 
catalogued emissivities, which were again based on solar 
photospheric abundances. The ratios between the pre¬ 
dicted and the measured fluxes provided the abundances 
relative to Fe, A/A(Fe), with respect to the correspond¬ 
ing solar ratios. In some cases, we measured the fluxes of 
more than one line for a single element. In these cases, 
we calculated the abundance for each line and computed 
averages, using \/cf\ as weights , where a a is the error in 
the abundance (see Sect. 15.3.31) . We implicitly assumed 
here that the abundance of an element is the same at 
all temperatures, an assumption th at is not necessaril y 
supported from solar observations l|, Iordan et alJll99'^ . 
The data quality at hand does not allow for further dis¬ 
crimination, however. 

Finally, we once more iterated the adjustment of the 
continuum level to obtain the absolute Fe abundance, as 
described above, now using the abundances determined 
from our procedure. The updated values closely agreed 
with the previously obtained Fe abundances. The error 
of the absolute Fe abundance was derived by varying it 
around its best ht value and requiring that the continuum 
fit be acceptable within one sigma. The final abundance 
values we report in this pap er refer to the solar photo- 
spheric abundances given bv I Anders k, Grevessel 1119891) . 
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except for Fe for which we adopt the value given by 
iGrevesse fc Sauvall l|l999H. 

5.3.3. Errors 

Errors arise from different sources. First, they are due 
to uncertainties in the atomic databases; these are not 
easily quantifiable, but are likely to be in the range of 
several percent to perhaps 20%, depending on the line 
under consideration. For method 2, we chose mostly 
bright, well-studied lines, and for the sake of definition we 
have assumed systematic uncertainties of 10% for each 
line. Statistical errors also arise from the fit of the 6- 
line model used to extract the line flux; these essentially 
originate from photon statistics. For each line flux, these 
two errors were summed in quadrature, and we call them 
line-errors a. 

The errors in the EMD were estimated by statistically 
varying the fluxes of the line blends according to their a, 
and repeating the EMD reconstruction for 19 different, 
perturbed line-flrpc lists. We thus derived formal upper 
and lower 1-sigma ranges of the EMD solutions by using 
standard formulae from Gaussian statistics, although we 
mention that the various solutions are not necessarily 
normally-distributed in log EM at any given tempera¬ 
ture. The standard deviations thus derived, however, 
provide a well-defined characteristic width of the distri¬ 
bution, and we verified that the ranges containing 68% 
of the solutions and the 1-sigma ranges are very simi¬ 
lar. Because the EM values scatter considerably in any 
temperature bin, we performed the statistics using log¬ 
arithms of the EMs in each temperature bin in order to 
avoid the average being biased by one or a few large val¬ 
ues. We extended the error analysis to include up to 
100 perturbed line lists but the error ranges did not sig¬ 
nificantly change. We caution that the EMD slopes on 
both sides of the peak temperature could be slightly de¬ 
pendent on the initial EMD guess, where we assume a 
slope of ±2 (see Sect. 15. 3. Til . We study the EMD results 
starting from different initial conditions for the slopes in 
Sect. 17.31 The final EMD slopes converged to similar val¬ 
ues. These effects are not taken into account in the EMD 
errors reported here. We note that these errors are only 
given as an indicator for the uncertainty in the EMD, but 
they are not explicitly used further in our error analysis. 

The error of the abundance A is proportional to the 
line-error, a a = A ■ a(x) / F(x). If the abundance is a 
weighted average, then the final error is the larger of i) 
l/cr^)“^/^ and ii) the error of the weighted means of 
all abundance values. Moreover, there is an error in the 
abundances arising from the variation in the 20 different 
EMD reconstructions: for each of these reconstructions, 
we found slightly different abundances with new errors. 
We defined the error from this variation as the larger of 
i) the average error found in each reconstruction, and ii) 
the standard deviation of the twenty abundance values 
per element. We note, however, that we adopted the 
abundance derived from the best-fit solution, with no 
perturbation applied. 

Finally, the error in the adjustment of the continuum 
(required to determine the line fluxes, Sect. 15. 3. 211 also 
affects the abundance errors. As our final error for a 
given abundance, we summed in quadrature the error 
related to the line-error (or the average if multiple lines 
were used, as defined above), the error arising from the 


variation of the EMD, and the error from the continuum 
adjustment. We emphasize, however, that this proce¬ 
dure can provide no more than a simulated estimate of 
realistic errors. The unknown systematic deviations in 
the atomic physics parameters prevent us from obtain¬ 
ing better estimates. 

6. RESULTS 

6.1. Emission Measure Distributions 

The EMDs derived from the two different methods 
are shown in Figure |S1 In the left column, EMDs from 
method 1 using SPEX in combination with a fit based 
on Chebychev polynomials of degree 6 and, where pos¬ 
sible, 8 are plotted. The middle and the right columns 
show EMDs reconstructed with our method 2, based on 
MEKAL and APEC emissivities, respectively. In the 
middle and right columns, the black histograms illustrate 
the best-fit EMDs while the red histograms mark the Icr 
range at each temperature, derived from the perturbed 
flux lists. As the best-fit EMDs are derived from unper¬ 
turbed fluxes, and they are not equal to the mean EMDs 
derived from the perturbed flux lists, these ranges of vari¬ 
ation do not need to be symmetrically arranged around 
the best-fit solutions. In some cases (EK Dra and 47 
Cas B), the lower error ranges drop rapidly to very low 
values at certain temperatures. Although this is a con¬ 
sequence of the increasing uncertainty in the EMD at 
the lowest and the highest temperatures, we note that 
the error ranges are given on a logarithmic scale; once 
the ranges becomes large, the precise level of the lower 
bound is of little importance. 

For method 2 the quality of the EMD can be mea¬ 
sured by comparing the predicted and the observed line 
fluxes. The final agreement between predicted and ob¬ 
served line fluxes is illustrated in Fig. El where we show 
the fractional deviation of the predicted line fluxes from 
the observed values, {Fq — F)/F for Fe and O Viii/0 vii 
flux ratio. Most line fluxes agree within 10-20%, with the 
larger deviations mainly relating to the weakest lines, i.e., 
the lines formed at high temperatures in the least active 
stars (e.g., Fexx for tt^ UMa). 

The EMDs derived from the different methods show 
rather similar characteristics. We see that the tempera¬ 
ture where the EMD peaks decreases toward older, less 
active stars, namely from about 10 MK for 47 Cas B and 
EK Dra to 5 MK for UMa, Ori ^tnd Cet, and to 
<4 MK for our oldest target, (3 Com. Characteristic val¬ 
ues for the Sun are 1-3 MK depe nding on the phase of its 
activity cycle llPeres et al.ll200(Tll . The average tempera¬ 
tures derived from method 2, logT, are listed in Table 13 
To obtain these values, we calculated the mean of logT, 
using the EMs in each bin as weights. Also given are the 
lower and upper threshold temperatures that comprise 
90% of the total EM (on each side of the EMD peak). 

With method 2, we generally obtain a smoother 
and flatter EMD than with method 1 and the Cheby¬ 
chev polynomial approximation. At temperatures above 
logT Ri 7.5, the EMD is not well constrained. This is 
obvious for the cooler coronae which do not provide any 
useful spectral lines at those temperatures, and we there¬ 
fore did not extend our DEM analysis to this range. Con¬ 
siderable scatter is still also found for the EMDs of 47 
Cas B and EK Dra, despite the availability of Fe xxiii- 
Fexxv lines. The reason resides in the fact that Fexxiii 
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and Fexxiv show very faint lines, and Fexxv is the 
only blend complex that covers the temperatures above 
logT Ri 7.5. In the more active stars, the EMDs from 
method 1 seem to be composed of two peaks at about 6 
and 20 MK, well separated by a local minimum. Also, 
in most EMDs derived from method 1, we hnd a deep 
decrease in the EM below about logT « 6.2, combined 
with a local EM peak around logT = 6.0. Both features 
become stronger if higher polynomial degrees are used. 

Several effects may contribute to this: First, a large 
range of solutions may in fact be compatible with 
the spectra, given that the spectral inversion is ill- 
conditioned, in particular in temperature regions where 
few constraints are available. Second, the fit of method 1 
iterates to minimum which considers only the Poisso- 
nian errors in the count spectrum and which may intro¬ 
duce EMD features of little relevance given the system¬ 
atic uncertainties in the atomic physics, while method 2 
has been terminated according to (see Eq. EJ, which 
approximately considers the atomic physics uncertainties 
as well. Third, method 1 uses many lines that may in¬ 
troduce uncertainty to the spectral ht, while the result 
of method 2 almost uniquely relies on the Oviii/Ovii 
flux ratio for the coolest portion of the EMD. And h- 
nally, method 1 imposes polynomial constraints on the 
solution, which favours the appearance of peaks and val¬ 
leys in the EMD, while method 2 starts with a smooth 
EMD that is changed only in so far as the spectrum re¬ 
quires. Eor example, if a line requires excess EM due to 
an underestimation of its emissivity at a given temper¬ 
ature, then the reconstruction process may compensate 
by lowering the EM at adjacent temperatures as dictated 
by lines dominating there. To test this hypothesis, we 
iterated method 2 excessively, to = 0.5. As the two 
examples in Eig.^lshow, very similar features also evolve 
in these examples. The low-temperature slopes appear 
to become partly steeper as well. 

We also note that the amplitudes of the oscillations 
are compatible with the error ranges from perturbing 
the line-flux lists used in method 2 (cf. Eig. O. It is 
conceivable that the oscillations found after a deep iter¬ 
ation of method 2 correspond to those seen in method 
1 although this cannot be explicitly proven, given the 
largely different approaches. However, the magnitudes 
of the oscillations appear to be similar. It is also possi¬ 
ble that the oscillations are present in the stellar EMD; 
we cannot reliably discriminate between this hypothesis 
and a numerical effect as long as we include statistical er¬ 
rors and assume the presence of systematic uncertainties 
of the magnitude adopted here (see Sect. 15.6.611 . 

In contrast, the abundance ratios turn out to be ro¬ 
bust, with no significant change when deeper iterations 
are applied. In fact, the synthesized spectra for the two 
cases are very similar, i.e., the two EMDs represent the 
spectra almost equally well. Comparing the synthesized 
spectrum with the observations in the wavelength inter¬ 
vals illustrated in Eig.|3 we find, for 47 Cas B, a reduced 
X^, Xred = 1-28 (for 1091 d.o.f) for the deeper integra¬ 
tion and Xred = 1-87 (for 1091 d.o.f) for our standard 
convergence criterion. Although at first sight this differ¬ 
ence appears significant, the important line features and 
the continuum in the RGS are htted well in both cases. 

We note that the EPIC MOS portion of the synthetic 


spectrum from method 2 is not very well fitted: the fit 
lies systematically below the data. This feature is prob¬ 
ably to be ascribed to a cross-calibration inaccuracy in 
the effective areas of the RGS and the EPIC MOS in¬ 
struments. As a matter of fact, a related effect is present 
in the spectrum from method 1 as well, but there, the 
continuum level of the synthetic spectrum is slightly but 
systematically too high in the RGS compared to the data, 
while the fit in the MOS spectrum is better. This is one 
important source for the somewhat higher Xred fo'^ 
method 2. We note, however, that the analysis based on 
method 2 corrects for the continuum discrepancy before 
line-flux extraction by adjusting the continuum level in¬ 
dividually both for the RGS and the EPIC spectra (see 
Sect. IH~6.III . 

In Table IHl we compare the Xred values with respect to 
the observations, for the synthetic spectra constructed 
from the EMDs that were obtained with method I and 
with method 2, respectively. Again, only the regions 
listed in Table 0 are used. The degrees of freedom are 
also listed in Table |S1 Although formally the same wave¬ 
length intervals were used, the number of degrees of free¬ 
dom are somewhat different in SPEX and XSPEC. This 
discrepancy comes from two different sources: first, par¬ 
tial bins at the beginning and the end of each interval are 
considered differently (XSPEC ignoring all partial bins). 
Although this discrepancy could be reduced by minor 
adjustments, we prefer to keep with the simple prescrip¬ 
tions of Table0]for easy reproduction of our results. The 
partial bins occur in relatively low-flux, shallow regions 
of the spectrum where essentially continuum is fitted, 
hence a difference by a single, usually well fitted con¬ 
tinuum bin is of little relevance. Second, the standard 
spectral software packages treat grouped bins that may 
contain bad data channels differently (SPEX breaks bins 
up into partial bins, while XSPEC does not). This is a 
feature of the standard software packages that we test 
here. 

The Xred from method 2 are slightly larger but the 
differences are not very substantial despite the system¬ 
atic uncertainties in the emissivities adopted in method 
2 but not in method 1. Beside the fact that in the re¬ 
constructed spectrum of method I the MOS spectrum 
is better fitted (given its higher S/N ratio), the slightly 
better Xred fo^ method 1 could also be a result of exces¬ 
sively deep iterations that aim at fitting poorly-described 
line fluxes at the cost of smoothness in the EMD. It is 
noteworthy that the Xred values from the deep iterations 
in method 2 closely approach the Xred values of method 
1 . 

There is ample literature on EMDs of active stars 
available. Several other authors have also found EMDs 
with features resembling ours. In particular, double- 
peaked EMDs ha ve previousl y been reported, e.g., by 
iMewe et al.llll996l) for AB Dord Ka astra et all lll996all for 
RS CVn binaries. iGiidel e t alJ lll997all for solar analogs, 
lSanz-For ca(^^ ^1; ^1. | l|200 1ll fo r the RS CVn binary A 
And. and iHuenemoerder et al.l ll2006ll for AR Lac, using 
entirely different EMD reconstruction methods. In the 
light of our discussion on errors and iteration depths in 
Sect.Elandini we cannot be certain on the actual reality of 
bi-modal DEMs, and given that other authors use similar 
atomic physics databases, the same caveat may apply to 














Coronal Evolution of the Sun in Time 


11 


other EMD reconstructions as well. iGiidel et al.l lll997all 
argued, based on solar-flare data, that a bi-modal struc¬ 
ture can arise from the rapid decay of the EM of a pop¬ 
ulation of flares as they are cooling. The present quality 
of our EMD inversions does not allow us to make more 
definitive conclusions at this point. 

6.2. Abundances 

The abundances found with the different methods are 
listed in Table 121 and plotted as a function of the FIP 
in Figure IHl We plot the abundance ratios A/¥e with 
respect to the solar photospheric ratios as a function of 
the FIP. The open circles represent the coronal abun¬ 
dances derived from method 1 whereas the filled circles 
show the abundances derived from method 2. We find 
a good agreement between the abundance sets, and an 
acceptable overall agreement between the results using 
the MEKAL database and the APEC database. 

Nevertheless, some differences can be noted in Table 121 
In the older stars Ori, tA Cet, and (3 Com, some 
systematic differences occur for the C/Fe abundance ra¬ 
tios, the abundances derived with APEC being smaller 
than those derived with SPEX although the error ranges 
are large. For the same stars, some differences are also 
present in the N/Fe abundance ratio. In both cases, the 
abundances are determined from only one faint line (C vi 
and N vii, respectively), making the measurements of the 
line fluxes difficult (see Tables |Sl and O . Differences also 
occur in the four older stars for the Ne/Fe abundance 
ratio. This is partly due to the lack of strong and re¬ 
liable (i.e., unblended) Ne lines. While in hot coronae, 
the strong Nex line serves as a reliable indicator for the 
Ne abundance, this line is much weaker in cooler coro¬ 
nae and strongly blended with Fe lines. The Neix lines 
at 13.55-13.7 A are always blended with Fe lines, but 
also become quite faint in the less active stars. Fur¬ 
ther, the inferred Ne flux in the line feature depends on 
the absolute Fe abundance. The latter is poorly deter¬ 
mined in particular in the less active stars where almost 
no continuum is present. We emphasize that we cannot 
attribute the discrepancy to any of the methods. The 
available data quality simply makes the determination 
of the Ne/Fe abundance ratio in cooler coronae ambigu¬ 
ous. 

Fe blending leads to some differences in the Mg/Fe 
abundances as well (the largest deviations are for tt^ UMa 
using APEC method 2, and Cet using SPEX method 
2). Finally the weakness of the S lines in the more active 
stars leads to some differences in the S/Fe abundances. 

The absolute Fe abundance was systematically higher 
when the APEC database was used. Note that in the 
spectrum of (3 Com, the continuum is almost nonexistent 
and the derivation of the Fe abundance is difficult. For 
the latter target, the d-fit did not converge for the Ne 
line with APEC. For this reason, these two points are 
missing in Tableland Figure El 

However, we recognize many of these features to be due 
to limitations of the data and the reconstruction meth¬ 
ods. On the other hand, these systematic differences are 
small compared to the general trends. As shown in Fig¬ 
ure El the abundances resulting from the two methods 
and the two databases agree mostly quite well within the 
errors, and the general trends are the same, regardless of 
the method used. 


6.3. Light curves 

The light curves of the six stars are shown in Figure El 
For each star, four light curves are plotted. They de¬ 
scribe, from top to bottom, the total count rates in the 
0.2 — 10 keV range (black), in the soft band (0.2 — 1 keV, 
green), in the hard band (> 1 keV, third plot), and the 
ratio between the hard and the soft count rates (blue). 
The upper energy thresholds used for the hard band vary 
from star to star and are listed in Table El Only data 
from detectors that operated in imaging mode were con¬ 
sidered (for 47 Cas B, EK Dra and (3 Com, the data from 
MOSl, MOS2 and PN were used; for tt^ UMa, data from 
MOSl and MOS2, and for Ori and Cet, only data 
from one MOS camera, MOS2 and MOSl respectively, 
were used). 

The light curves display considerable variability. We 
observe the presence of large flares on all stars in our sam¬ 
ple except (3 Com. However, even after excluding these 
flares, the light curves still show considerable variability 
that cannot be described by steady, quiescent coronal 
emission. We also note that the hard emission becomes 
weaker toward older stars, in agreement with the decline 
of the average coronal temperature described earlier. 

7. DISCUSSION 

7.1. Correlation between the Parameters 

The results of our analysis of six solar analogs clearly 
show a number of trends that we wish to quantify be¬ 
low using our results from method 2. Before doing so, 
we note that the X-ray results from the observation of 
Cet are rather similar to results from tt^ UMa and Ori, 
despite the former’s significantly longer rotation period 
and, hence, higher age. The reason for this discrepancy is 
not entirely clear, but we note that Cet has a slightly 
later spectral type and therefore a somewhat lower mass 
than the other tar gets. It is known that l ater-type stars 
evolve more slowly llSoderhlnm et alJl993fl and that they 
remain in a state of maximum X-ray lum inosity (the sat¬ 
urati on limit) for longer rotation periods llPizzolato et alj 
Hoi than stars of earlier spectral type. Both effects 
make Cet look somewhat younger than inferred from 
a rotation-activity relation that is appropriate for early 
G stars. 

_ In Figure EDI we plot the mean coronal temperature, 
T, (see TableEJ as a function of the total luminosity Lx- 
We fitted the data with a power-law. Because both vari¬ 
ables to be correlated are likely to be affected by system¬ 
atic scatter around any power-law, we use t he ordinary 
least s quares bisector method as described bv llsobe et alJ 
iTT^ . Clearly, the two parameters are correlated. We 
find, for our results from MEKAL and APEC, respec¬ 
tively. 

Lx « 1.17 X 1026t4-26±o.4i ergs“i (MEKAL), (6) 

Lx«1.61 X ergs"^ (APEC), (7) 

where T is i n MK. These relatio ns are consistent with 
the results of iCiidel et alJ l|1997ar . except that the latter 
authors used the higher temperature for a model with 
two thermal components fitted to ROSAT spectra. Two 
points for the Sun at minimum and maxim um activity 
level are also plotted for comparison (after iPeres et alJ 
|2003). 
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In Figure CH the temperature as a function of the pe¬ 
riod is shown. Again, we fitted the data with a power- 
law. In this case, the relations calculated from the two 
power-laws are given by 

MK (MEKAL), (8) 

MK (APEC), (9) 

where Pmt is the rotation period in days. Similar re- 
sults were obtained for the ROSAT data bv lGiidel et al.1 
dnizi, again considering the higher fit-temperature in¬ 
stead of the mean temperature. As we averaged the tem¬ 
perature with the EM in each bin used as weights, the 
steady decrease of the mean temperature with period is 
consistent with a decrease of the amount of hot plasma 
as the star spins down. 

Equations mni allow us to check for consistency with 
published relations between Lx and Prot- By combining 
formula (6) with (8) or, respectively, (7) with (9), we 
find 

Lx = 4.01 X erg s’^ (MEKAL) j(10) 

Lx = 4.04 X 103°P-2 °3±0.35 erg s-i (APEC), (11) 

which is cons istent with previously reported dependence s 
of this type llPallavicini et al.lll98lHCudel et alJll997all . 
A linear regression for log Lx and log Prot yields the same 
result, with a power-law index of -2.03. 

We also studied the relation between radio luminosity 
and the temperature. Radio luminosities or upper limits 
thereof are available l|Ciidel fc Caidoijl200lH for five out 
of the six targets and are plotted in FigurelT^ They refer 
to low emission levels outside obvious flares. We use the 
values of the upper limits in the regression analysis. The 
slopes of the power-laws are therefore lower limits, given 
by 

Lr R:! 0.86 X 109t 3.65±0.46 g-1 jj^-l (MEKA^)?) 

LfiWl.69 X erg s"^ Hz-1 (APEC)(13) 

These relations suggest a relation between the nonther- 
mal electron population, responsible for radio gyrosyn- 
chrotron emission, and coronal heating. 

7.2. Abundances 

In the solar corona, the so-called FIP effect has been 
observed, in which the elements with a FIP lower than 
10 eV are overabundant relative to the solar photo- 
spheric composition, whereas the elements with a higher 
FIP show the same abundance as the solar photosphere 
^el dma.nlll 992t iLa.ming et aUlTo^ iFeldma.n L: Laminel 

1200011 . Recent spectroscopic analysis with XMM-Newton 

and Chandra has shown that in very active stars, an 
inverse effect is present, in which the low-FIP ele- 
ments are depleted re lative to the high-FIP elements 
llBrinkman et al.l 1200111 . In our sample, we observe an 
evolutionary trend from an inverse FIP effect for the most 
active star 47 Cas B to a solar-like FIP effect in the old¬ 
est stars (FigurelHl). We note, however, that the absolute 
abundances of low-FIP elements such as Fe do not reach 
values as high as in the solar corona, where overabun - 
dances by factors of a few are common llFeldmanlll992^ . 

In Figure ESI the abundances of Fe, Ne, and the 
ratios A(Ne)/A(Fe), A(0)/A(Ne), A(0)/A(Fe), and 
A(Mg)/A(Fe) are plotted as a function of the temper¬ 
ature (based on MEKAL/SPEX, method 2). The dot¬ 
ted regions include the ranges of a larger stellar sample 


llCiidell 1200411 . The abundance of the low-FIP element 
Fe tends to decrease from a nearly photospheric value 
for stars with an average temperature of 3 to 5 MK to a 
lower abundance of « 0.5 for the two more active stars 
EK Dra and 47 Cas B. The error bars for the coolest 
star (/3 Com) are not plotted for Fe and Ne, since they 
exceed the range illustrated in the figures. This is due to 
the near-absence of a continuum in this star, which makes 
absolute abundance determinations difficu lt. The abun- 
dance ratios, however, are robust (see also lAudard et alJ 
1 ^ . 

In the middle left and the bottom left panels, the abun¬ 
dance ratio of Ne/Fe and 0/Fe, respectively, are shown 
as a function of the average coronal temperature. Be¬ 
cause Ne and O are high-FIP elements, and Fe is a low- 
FIP element, the ratios increase with temperature, con¬ 
firming the evolutionary trend of decreasing low-FIP el¬ 
ements (such as Fe) with increasing activ ity. A similar 
effect was found bv lAudard et al.l ^^2003a^ for a sample 
of active RS CVn-type binaries. This is also confirmed 
by a larger sample of active stars llCiidelll2nn4 . 

The plots in the middle right and bottom right pan¬ 
els show the abundance ratio of 0/Ne, two high-FIP ele¬ 
ments, and of Mg/Fe, two low-FIP elements. Our sample 
of stars is too small to constrain a trend i n these plots . 
However, the larger star sample studied bv lCiidell (|2004ll 
shows a nearly flat distribution for both ratios. 

Could it be that the coronal abundance pattern reflects 
the composition of the underlying photosphere? This 
view does not find support from other studies. Although 
a full picture would include knowledge of the abundances 
of other elements such as C, N, O, we see, from the sum¬ 
mary in Sect. o no support for a photospheric abun¬ 
dance pattern that significantly deviates from solar. Fi¬ 
nally, we recall that the solar coronal composition does 
not reflect the photospheric composition either, hence an 
agreement between photospheric and coronal abundances 
is not a priori anticipated for solar analogs. 

7.3. Flares and Coronal Heating 

In the previous sections, we have highlighted correla¬ 
tions between observable parameters, and we have found 
continuous variability in all six targets to an extent that 
hardly any time interval is free of fluctuations. Al¬ 
though conventional interpretation of coronal structure 
often makes use of the app roximation of static coronal 
loops (as, e.g., described bv iRosner et al ] lI978ll . the in¬ 
terpretation of the phenomenology revealed by our light 
curves cannot be rooted in strictly static loops, although 
static loop models may, under certain circumstances, 
serve as approximati ons even under flaring conditions 
Ijjakimiec et al.llI992ll. 

There is appeal in the alternative, extreme model as¬ 
suming that the coronal emission is entirely due to dy¬ 
namic, flaring loops. A number of observed features re¬ 
ported in this paper and in the previous literature seem 
to support such a model: i) with the sensitivity avail¬ 
able with XMM-Newton and Chandra, X-ray emission 
previously ascribed to a quiescent component is now rec¬ 
ognized to be continuously variable; in the most extreme 
cases, no steady co mponent can reasonab ly be identified 
in the light curves l|Audard et al.l[^03b|l . ii) More ac¬ 
tive stars (i.e., stars with a higher Lx/Lboi ratio) appear 
to maintain hotter coronae. This is difficult to explain 
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with a model that assumes a corona composed only of 
steady magnetic loops, with the principal determinant of 
Lx being the magnetic filling factor, up to the empirical 
saturation value of logLx/^boi ~ —3. Such models do 
not automatically explain why the coronal temperature 
increases with increasing Lx- We discuss this point fur¬ 
ther in our conclusions below, iii) Active stars continu¬ 
ously produce radio emission from accelerated electrons. 
The lifetime of the latter is short, probably amountin g 
to no more than seconds to minutes llKundu et al Inna). 
In the solar corona, flare energy-release processes are re¬ 
quired to produce such electron populations. In order to 
generate the observed stellar radio emission, high-energy 
electrons must be replenished frequently. 

This question has previously been addressed by study¬ 
ing, for a model in which the corona i s heated en¬ 
tirely by flares, the EMD llGiidell lIQQTd: iGiidel et al.1 
12003 ), th e average temperatu re to be ascribed to such 
a corona l|Audard et al.ll200nll . and the light curve char- 
acteris tics expected from a supe r position of stochastic 
flares llAudard et al.l 119991 120001: jlOagl waD et al ] i2ooa 

IGiidel et al.1120031: lArzner fc Giidell 1200411 . We are now 

in a position to apply the methodology to our results. 

First, we use the shapes of our EMDs to character¬ 
ize the underl ying flare populatio n in the framework 
of this model. IGiidel et al.l l|2003ll derived an analytic 
expression for a DEM of a flare-heated corona under 
the assumption that the temperature and the flare den¬ 
sity both decay exponentially with time constants tt 
and Tn, respectively, and that a relation between flare 
peak-EM and peak-T holds, EMp oc T^, as reported by 
iFeldma.n et alJ ill Dfllifl . Then, the DEM follows power- 
law relations on both sides of its peak temperature (Tm), 
with 


VG 1 O'- I rp-(b-(j>){a-2f3)/{l-0)+2b-(j> rp ^ rp^ 

™(14) 

where b « 4.3 ± 0.35 as de rived from a large sample 
of stellar flares llGiidelll2nn4(l . The temperature Tm de¬ 
pends on the energy of the smallest flares participatin g 
in the heating in this simple model (IGiidel et al.ll200^ . 
Further, /3 is a power-law index for a relation between 
the flare e-folding de cay time and its r adiated energy, 
T (X . As noted bv IGiidel et al.l ll2003ll . (3 is probably 
close to zero although an extreme case oi j3 = 0.25 was 
also studied. The variable (j) gives the slope of the cool¬ 
ing function (radiative power per unit emission measure) 
for a power-law approximation in the temperature inter¬ 
val of interest. We use (j) = —0.3 as an approximation 
of the slope of the cooling function in the logarithmic 
temperature interval 6.8-7.5, i.e. the temperature in- 
terval above the EM D peak temperature (see Fig. 10 in 
lAudard et alJl2C)(i^ . The parameter a is of primary in¬ 
terest for us: It is the exponent of the distribution of 
the occurrence rate N of flares in radiated energy, viz., 
dN/d E oc E~°‘ as found for a large sample of solar flares 
(e.g., iCrosbv et al.lll99^ . but applicable also to stellar 
coronae (see references above). In the first equation, ap¬ 
plicable to the cooler portion of the DEM, C = Tn/rT- 
This parameter describes the amount of heating occur¬ 
ring during the flare decay, with C = 2 corresponding to 
free cooling without heating, and C ~ 0.5 co rresponding 
to extreme heating rates during the decay llReale et al.l 


[1993). 

We have measured the slopes of our EMDs on both 
sides of Tm- To find the possible ranges for the best-fit 
slopes, we re-analyzed our data using method 2 by start¬ 
ing with different initial conditions for the slopes in our 
iteration, but the DEMs converged to similar values. The 
ranges of the resulting best fit slopes (not considering the 
error ranges in the DEM) are reported in Table 11(11 We 
then used Eq. nn to determine the most likely C, and a 
with their acceptable ranges. The results are also re¬ 
ported in Table nni 

From the low-T slopes of the EMDs derived with 
method 2 using SPEX, we find ^ to be around unity 
in all cases. Such values are typical for individu al flares 
observed on active stars (see IGiidel et al.11200.'^ and ref¬ 
erences therein) and support our assumption that such 
flares contribute significantly to the overall observed 
emission. Static loops, on the other hand, normally 
produce shallower DEMs, with slopes of -1-1 to -1-3/2, 
depending on t he amount of conductive flux at the 
loop footpoints llRosner et alJll97^ Ivan den Oord et alJ 

|T993- 

From the high-T slope, we derive a « 2.2—2.8, in excel¬ 
lent agreement with a values determined from long light 
curves of active stars llAud ard et alJ200ftl KashvaD et alJ 

l2flfl2tlGiide1 et al Jl2nn.3n^^ue^! Giidelll2nn4|) . If rr > 2 

and the emitted energy is integrated over the flare- 
rate distribution to obtain the total radiative loss, i.e., 

-btot = E{dN/dE)dE, then Ltot diverges as Ei 0, 
i.e., the smallest flares dominate coronal heating, and a 
lower energy cutoff is required for this power-law. 

We now continue this consideration by simulating light 
curves based on a characteristic flare shape for the 
three most active stars, assuming a flare-rate distribu¬ 
tion based on a as determined above. The flare shape 
was derived from the convolution of an exponential func¬ 
tion (cut off at t < 0, describing the decay) and a Gaus¬ 
sian (important to describe the rise and the peak phase). 
The shapes of the largest flares in the light curves were 
used to estimate the characteristic rise and the decay 
time parameters. We performed two sets of simulations: 
one with (3 = 0 and one with /3 = 0.25. In each simulation 
a was chosen within the acceptable range for a given (3. 
The flare decay-time was varied according to the scaling 
T cx E^. For the largest flare in each light curve, we mea¬ 
sured the amplitude and fitted the decay phase with an 
exponential function to find the decay time, and hence 
the emitted energy. We then set the maximum flare en¬ 
ergy E 2 equal to the energy of the largest observed flare 
and constructed a power-law distribution of flares in en¬ 
ergy, down to a selectable minimum flare energy Ei. The 
flares were randomly distributed in time (assuming a to¬ 
tal of 40 ks), and their light curves were superimposed. 
The rate of flares at a given energy, equivalent to the 
probability of the largest flares to occur within the sim¬ 
ulation time, could be statistically varied. We measured 
the modulation depth, i.e., the ratio between the root- 
mean-square scatter of the light curve and the average 
luminosity, and compared it with the same measure for 
the observed light curves outside the outstanding largest 
one or two flares. The minimum flare energy Ei and 
the numbers of large flares occurring was varied until 
the modulation depth and the average luminosity level 
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agreed with the observed light curves. The modulation 
depth of the observed light curves is given in Table El 
It increases with decreasing activity. The simulated light 
curves, for /3 = 0 and a = 2.25, 2.28, and 2.54 for 47 
Cas, EK Dra, and UMa respectively, are shown in 
FigureEl The range logi? 2 /ifi (in dex) of flare energies 
thus required for the light curve is also given in Table 11171 

We see that typically 1.8—5.3 orders of magnitude of 
flare energies for /3 = 0 and 1.8—6.8 for/3 = 0.25, respec¬ 
tively, are required to “describe” our light curves. The 
interesting point is that this range is larger for the more 
luminous stars. This can be understood because a nar¬ 
row range of flare energies produces a more strongly mod¬ 
ulated light curve, while a large number of small flares 
merely adds a quasi-steady baseline level. However, in 
Equation El the turnover in the DEM is determined by 
the smallest flares participating in the statistics. One 
would thus expect that Ex decreases with decreasing ac¬ 
tivity, contrary to the results in Table El However, we 
emphasize that our light curve model is based on the ex¬ 
treme assumption that all emission originates from flares. 
If there is a steady baseline level of X-ray emission not 
directly related to flare decays, as suggested from solar 
observations, we can obviously not determine Ei from 
light curve analysis. Also, our model does not take into 
account any possible change in a at lower energies. We 
know from solar observations that there will be a contin¬ 
uation toward smaller flares as well, rather than a lower 
energy threshold. The values in Table El are therefore 
only indicative of the range of flare energies required in 
the most extreme model discussed here. 

8. CONCLUSIONS 

As the rotation rate of a solar analog decreases dur¬ 
ing its evolution on the main sequence, the efficiency of 
the internal dynamo weakens, resulting in a decrease of 
the magnetic activity in the stellar atmosphere. We have 
studied systematic trends in the long-term evolution of 
stellar coronal X-ray emission for ages in the range 0.1- 
2 Gyr. As a consequence of the stellar spin-down, the 
X-ray luminosity steadily decreases from levels that may 
be close to the empirical saturation limit in the youngest 
stars (Tx/Tboi ~ 10“^), to levels approximately two or¬ 
ders of magnitude lower within the hrst two Gyr. During 
the next w3-4 Gyr, Lx reduces by another factor of ten 
to leve ls as seen in a Gen, f3 Hyi, or the Sun llGiidel et al.1 
noOTall . The overall trends reported here conhrm ear- 
lier studies based on lower-resolution X- ray spectroscopy 
llMaggio et al.lll98'^ iGiidel et al.lll997^ . 

The high-resolution X-ray spectroscopy now available 
has permitted a more detailed study of the composition 
and the thermal structure of solar-like coronae than was 
previously possible. We have studied emission measure 
distributions and coronal element abundances for all six 
targets in a homogeneous way, applying two widely dif¬ 
fering methods and using two different sets of atomic 
parameters. It is important to recall that both of our 
methods include modeling of line blends to the extent 
possible with the presently available line emissivities. 

There is gratifying agreement between the results from 
the two methods, although it appears that the choice of 
the atomic database introduces systematic differences. 
This is perhaps not surprising as the presently avail¬ 
able compilations of atomic parameters are incomplete 


and suffer from systematic uncertainties. A more serious 
limitation is set by the mathematical problem of spec¬ 
tral inversion itself. While even counting statistics at the 
percent level makes th e inversion pro blem ill-conditioned 
ilGraig fc BrownI [197(11 1, fudge 1 1200211 . systematic uncer¬ 
tainties in the line emissivities of perhaps up to 10-20% 
may introduce various structure in our EMDs that may 
not correspond to coronal features. We have suggested 
that a reasonable convergence criterion should be set, 
although it is difficult to assess at what level artihcial 
structure is introduced into a reconstructed EMD. Re¬ 
gardless of these inconsistencies, however, we recover el¬ 
ement abundances that are rather robust. 

With these limitations in mind, we have found system¬ 
atic trends in the EMD structure as a star ages. Not only 
does the total emission measure continuously decrease, 
the temperature where the EMD peaks also decays with 
time. The EM-weighted logarithmic average of the coro¬ 
nal temperature T thus follows a power-law dependence 
on the X-ray luminosity, namely T oc A 

similar relation between the dominant c oronal temper¬ 
ature and Lx was already studied by iSchriiver et alJ 
jUlL using Einste in data of a large sample of stars. 
iSchmitt et al.l iHQQflU found f or l-T model fits T (x 
from Einstein data. Einally, IGiidel et a 0 111 ffi)7all used 
ROSAT data and a G-star sample similar to ours and 
found a relation between the higher temperature of a 2- 
T model and Lx in complete agreement with our results. 

The observed trends cannot be explained by a model 
that is based exclusively on different filling factors of 
the surface magnetic field, as an increased filling factor 
does not explain per se why the temperature should in¬ 
crease. We have discussed an extreme case of an alterna¬ 
tive model in which a statistical distribution of flares is 
responsible for the correlation between Lx and T. The 
question then is why the flare rate (above any given base- 
level energy) is higher in more active (i.e., more rapidly 
rotating) stars. Our data cannot give a conclusive an¬ 
swer. One possibility is that more active stars indeed 
do show a larger magnetic surface filling factor, and that 
the higher density of magnetic loops leads to more mag- 
netic reconnection, thus producing a higher flare rate 
ijGiidel et al.l 11997^ : this, then, also includes a more 
prominent population of very large flares that produce 
both large emission measures and very high tempera¬ 
tures, thus shifting the average temperature to higher 
values. Once the magnetic filling dilutes, the interac¬ 
tions between neighboring magnetic loop systems will 
become less freque nt, and both Lx and T decrease (see 
IGiidel et aLlll997all. 

Alternatively, instead of increasing the surface filling 
factor, other ingredients may systematically change with 
changing rotation period, such as the structure of surface 
magnetic field in active regions or differences in the con¬ 
vection pattern that jostles the magnetic-loop footpoints. 
In a statistical-flare model, then, the way to explain the 
Lx — T correlation would be to reheat the same active 
regions more frequently in more active stars, as some pro¬ 
cess brings non-potential energy i nto the magnet i c field s 
at a higher rate. As argued by lAudard et al.l 1)200011 . 
such frequently-heated loops can be approximated by 
static loops although the heating process is non-static. 
In this case, approximating the loop temperature with its 
apex temperature (where most of the EM is found), the 
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iRosner et al loop scaling law predicts Lx cjc 

for loops of given length below the coronal pressure scale 
heig ht, and Lx (x for loops larger than this limit 
Isee lAudard et all 120001 for details). These predictions 
are close to our observational finding. On the other hand, 
static loops would predict slopes of the DEM on the low- 
temperature side that are significantly smaller than those 
determined by us. 

To conclude, we are presently unable to distinguish 
between a model in which the flare rate is controlled by 
the magnetic filling factor, and one in which constrained 
active regions flare progressively more frequently as the 
rotation rate of a star increases. Both approaches, how¬ 
ever, are compatible with the hypothesis that much of 
the coronal heating is induced by flaring, regardless of 
the ultimate cause of the increased flare rate in more 
active stars. 

Because lar ger _ flares produce hotter plasma 

llFeldman et alJ 1199511 . more active stars produce 
hotter coronae. This trend is unequivocally recovered 
from our observations and further supports a picture 
in which flares cont ribute significantly to the overall 
coronal heating (e g.. lAudard et alJl2000t IKashvao et al.1 
[2f)f)l ICiidel TtMlElf)^ 

We note in passing that an alternative view with 
similar con s equen ces has recently been presented by 
iPeres et "mI ll2004ll . These authors suggest that, based 
on properties of coronal structures seen on the Sun, a 
higher occurrence of very compact, hot features includ¬ 
ing flares make more active coronae hotter. 

We have also derived element abundances and found 
good agreement between our two methods. Somewhat 
more systematic deviations can be noted if different 
atomic databases are used (MEKAL, APEC), but the 
trends in the abundance pattern agree and markedly 
change with changing activity l evel. A similar tren d 
was noted for RS CVn binaries bv lAudard et all l|2003all . 
but the latter study referred to extremely active stars in 


Anders, E., &; Grevesse, N. 1989, Geochim. Gosmochim. Acta, 53, 
197 

Arnaud, K. A. 1996, in Astronomical Data Analysis Software and 
Systems V, ed. G. Jacoby &; J. Barnes (San Francisco: ASP), 17 
Arnaud, M. & Raymond, J. 1992, ApJ, 398, 394 
Arnaud, M. & Rothenflug, R. 1985, A&AS, 60, 425 
Arzner, K., & Giidel, M. 2004, ApJ, 602, 363 
Audard, M., Giidel, M., Drake, J. J., &; Kashyap, V. 2000, ApJ, 
541, 396 

Audard, M., Giidel, M., & Guinan, E. F. 1999, ApJ, 513, L53 
Audard, M., Giidel, M., & Skinner, S. L. 200313, ApJ, 589, 983 
Audard, M., Giidel, M., Sres, A., Raassen, A. J. J., & Mewe, R. 
2003a, A&A, 398, 1137 

Audard, M., Telleschi, A., Giidel, M., Skinner, S. L., Pallavicini, 
R., & Mitra-Kraev, U. 2004, ApJ, in press 
Brinkman, A. G., et al. 2001, A&A, 365, L324 
Bouvier, J. 1990, AJ, 99, 946 

Ganuto, V. M., Levine, J. S., Augustsson, T. R., & Imhoff, C. L. 
1982, Nature, 296, 816 

Gayrel de Strobel, G., Soubiran, C. & Ralite, N. 2001, A&A, 373, 
159 

Clegg, R. E. S., Tomkin, J., & Lambert, D. L. 1981, ApJ, 250, 262 
Craig, J. D., & Brown, J. C. 1976, A&A, 49, 239 
Crosby, N. B., Aschwanden, M. J., & Dennis, B. R. 1993, Solar 
Phys., 143, 275 

den Herder, J. W., et al. 2001, A&A, 365, L7 

Dorren, J. D., Giidel, M., & Guinan, E. F. 1995, ApJ, 448, 431 

Dorren, J. D., & Guinan, E. F. 1994, ApJ, 428, 805 


which the abundance pattern changed from a strong in¬ 
verse FIP effect to a flat distribution with decreasing ac¬ 
tivity (or mean coronal temperature). In our sample, the 
change from an inverse or flat distribution to a solar-like 
distribution occurs at ages of less than 300 Myr or rota¬ 
tion periods longer than « 3 days. Incidentally, a rapid 
decay of nonthermal radio emission has been noted for 
the same activity range. We hypothesize that the same 
electrons that are responsible for the observed gyrosyn- 
chrotron emission also induce an in verse-FIP effect in the 
most active stars, as follows (see iGiidel et al.l 120021 for 
further arguments): if electrons are streaming along the 
magnetic fields toward the chromosphere, they build up 
a downward-pointing electric field that acts to suppress 
positive currents from the chromosphere to the corona. 
In other words, ions in the chromosphere are prevented 
from streaming into the corona, while neutral, predomi¬ 
nantly high-FIP elements, are not affected. As the elec¬ 
tron population diminishes in less active stars, the sup¬ 
pression of ion diffusion into the corona disappears, and 
a solar-like FIP effect can build up . by wh a tever (still 
unidentified) mechanisms. Recently. iLamin^ ll2004f) pre¬ 
sented an alternative model in which both the solar-like 
FIP and the inverse FIP effect are related to a common 
plasma-physical cause. 

The authors acknowledge helpful comments by the ref¬ 
eree. This research is based on observations obtained 
with XMM-Newton, an ESA science mission with instru¬ 
ments and contributions directly funded by ESA Mem¬ 
ber States and the USA (NASA). AT and MG acknowl¬ 
edge support from the Swiss National Science Founda¬ 
tion (grant 20-66875) and from the Swiss Academy of 
Natural Sciences. MA acknowledges support from NASA 
to Columbia University for XMM-Newton mission sup¬ 
port and data analysis. SS acknowledges support from 
NASA/GSFC grant NAG5-13677. This research made 
use of the SIMBAD database, operated by CDS, Stras¬ 
bourg. 


Drake, J. J., Brickhouse, N. S., Kashyap, V., Laming, J. M., 
Huenemoerder, D. P., Smith, R., &; Wargelin, B. J. 2001, ApJ, 
548, LSI 

Drake, J. J., Laming, J. M., &; Widing, K. G. 1997, ApJ, 478, 403 
Drake, S. A., Singh, K. P., White, N. E., &; Simon, T. 1994, ApJ, 
436, L87 

Edvardsson, B., Andersen, J., Gustafsson, B., Lambert, D. L., 
Nissen, R E., Tomkin, J. 1993, A&A, 275, 101 
Feldman, U. 1992, Phys. Scripta, 46, 202 
Feldman, U. & Laming, J. M. 2000, Phys. Scripta, 61, 222 
Feldman, U., Laming, J. M., & Doschek, G. A. 1995, ApJ, 451, L79 
A&A, 377, 911 

Gagne M., Caillault J., & Stauffer J. R. 1995, ApJ, 450, 217 
Gaidos, E. J.,& Gonzalez, G. 2002, New Astron., 7, 211 
Gratton, R. G., Carretta, E., & Gastelli, F. 1996, A&A, 314, 191 
Grevesse, N., & Sauval, A. J. 1999, A&A, 347, 348 
Giidel, M. 1997, ApJ, 480, L121 
Giidel, M. 2004, A&AR, 12, 71 

Giidel, M., Audard, M., Kashyap, V. L., Drake, J. J., & Guinan, 

E. F. 2003, ApJ, 582, 423 

Giidel, M., Audard, M., Sres, A., Wehrli, R., Behar, E., Mewe, 
R., Raassen, A. J. J., & Magee, H. R. M. 2002, Stellar Coronae 
in the Chandra and XMM-Newton Era, Eds. F. Favata & J. J. 
Drake (San Francisco: ASP), 497 
Giidel, M., & Gaidos, E. J. 2001, in 11th Eleventh Cambridge 
Workshop on Cool Stars, Stellar Systems, and the Sun. Ed. 
R. J. Garcia Lopez, R. Rebolo, & M. R. Zapatero Osorio (San 
Francisco: ASP), 662 


























16 


Telleschi et al. 


Giidel, M., Guinan, E. F., Etzel, P. B., Mewe, R., Kaastra, J. S., 
& Skinner, S. L. 1998a, in Tenth Cambridge Workshop on Cool 
Stars, Stellar Systems and the Sun, eds. R. A. Donahue & J. A. 
Bookbinder (San Francisco: ASP), 1247 
Giidel, M., Guinan, E. F., Mewe, R., Kaastra, J. S., & Skinner, 
S. L. 1997b, ApJ, 479, 416 

Giidel, M., Guinan, E. F., & Skinner, S. L. 1997a, ApJ, 483, 947 
Giidel, M., Guinan, E. F., & Skinner, S. L. 1998b, in Tenth 
Cambridge Workshop on Cool Stars, Stellar Systems and the 
Sun, eds. R. A. Donahue & J. A. Bookbinder (San Francisco; 
ASP), 1041 

Giidel, M., Schmitt, J. H. M. M., & Benz, A. O. 1995, A&A, 293, 
L49 

Giidel, M., et al. 2001, A&A, 365, L336 

Guinan, E. F., Ribas, I., & Harper, G. M. 2003, ApJ, 594, 561 
Gray, R. O., Graham, P. W., & Hoyt, S. R. 2001, AJ, 121, 2159 
Huenemoerder, D. P., Canizares, C. R., Drake, J. J., Sanz-Focada, 
J. 2003, ApJ, 595, 1131 

Isobe, T., Feigelson, E. D., Akritas, M. G., & Babu, G. J. 1990, 
ApJ, 364, 104 

Jakimiec, J., Sylwester, B., Sylwester, J., Serio, S., Peres, G., & 
Reale, F. 1992, A&A, 253, 269 
Jansen, F., et al. 2001, A&A, 365, LI 

Jordan, C.,Doscheck, G. A., Drake, J. J., Glavin, A. B., & 
Raymond,J. C. 1998 in ASP Conf. Ser. 154, The Tenth 
Cambridge Workshop on Cool Stars, Stellar Systems and the 
Sun, Edited by R. A. Donahue and J. A. Bookbinder, 91 
Judge, P. 2002, Stellar Coronae in the Chandra and XMM-Newton 
Era, ASP Conference Proceedings, Vol. 277, ed. Favata & Drake, 
(San Francisco: Astronomical Society of the Pacific), 45 
Kaastra, J. S., Mewe, R., Liedahl, D. A., Singh, K. P., White, N. E., 
& Drake, S. A. 1996a, A&A, 314, 547 
Kaastra, J. S., Mewe, R., & Nieuwenhuijzen, H. 1996b, in UV and 
X-ray Spectroscopy of Astrophysical and Laboratory Plasmas, 
ed. K. Yamashita & T. Watanabe (Tokyo: Universal Academy), 
411 

Kashyap, V., Drake, J. J., Giidel, M., & Audard, M. 2002, ApJ, 
580, 1118 

Kundu, M. R., Jackson, P. D., White, S. M., & Melozzi, M. 1987, 
ApJ, 312, 822 

Laming, J. M. 2004, ApJ, 614, 1063 

Laming, J. M., Drake, J. J., & Widing, K. G. 1995, ApJ, 443, 416 
Laming, J. M., Drake, J. J., & Widing, K. G. 1996, ApJ, 462, 948 
MacGregor, K. B., & Brenner, M. 1991, ApJ, 376, 204 
Maggio, A., Sciortino, S., Vaiana, G. S., Majer, P., Bookbinder, 
J. A., Golub, L., Hamden, F. R. Jr., & Rosner, R. 1987, ApJ, 
315, 687 

Mazzotta, P., Mazzitelli, G., Colafrancesco, S. & Vittorio, N. 1998, 
A&AS 133, 403 

Mermilliod J.-C., Rosvick J. M., Duquennoy A., & Mayor M. 1992, 
A&A, 265, 513 

Mewe, R. 1999, in X-Ray Spectroscopy in Astrophysics, eds. J. van 
Paradijs & J. A. M. Bleeker (Berlin: Springer), 109 
Mewe, R., Gronenschild, E. H. B. M., & van den Oord, G. H. J. 
1985, A&AS, 62, 197 

Mewe, R., Kaastra, J. S., White, S. M., & Pallavicini, R. 1996, 
A&A, 315, 170 

Meyer, J.-P. 1985a, ApJS, 57, 151 
Meyer, J.-P. 1985b, ApJS, 57, 171 
Micela, G., & Marino, A. 2003, A&A, 404, 637 


Micela G., Sciortino S., Kashyap V., Hamden F. R., & Rosner R. 

1996, ApJS, 102, 75 
Micela G., et al. 1999, A&A, 341, 751 

Ness, J.-U., Giidel, M., Schmitt, J. H. M. M., Audard, M., & 
Telleschi, A. 2004, A&A, 427, 667 
Ness, J.-U., Schmitt, J. H. M. M., Audard, M., Giidel, M., & Mewe 
R. 2003, A&A, 407, 347 

Ottmann, R., Pfeiffer, M. J., & Gehren, T. 1998, A&A, 338, 661 
Pallavicini, R., Golub, L., Rosner, R., Vaiana, G. S., Ayres, T., & 
Linsky, J. L. 1981, ApJ, 248, 279 
Panzera, M. R., Tagliaferri, G., Pasinetti, L., & Antonello, E. 1999, 
A&A, 348, 161 

Paresce, F. 1984, AJ, 89, 1022 
Peres, G., Orlando, S., & Reale, F. ApJ, 612, 472 
Peres, G., Orlando, S., Reale, F., Rosner, R., & Hudson, H. 2000, 
ApJ, 528, 537 

Perryman, M. A. C., and the Hipparcos Science Team 1997, 
The Hipparcos and Tycho Catalogues (ESA report SP-1200; 
Noordwijk: ESA) 

Pinsonneault, M. H., Kawaler, S. D. , Sofia, S. & Demarque, P. 
1989, ApJ, 338, 424 

Pizzolato, N., Maggio, A., Micela, G., Sciortino, S., & Ventura, P. 
2003, A&A, 397, 147 

Raassen, A. J. J., Ness, J.-U., Mewe, R., van der Meer, R. L. J., 
Burwitz, V., & Kaastra, J. S. 2003, A&A, 400, 671 
Reale, F., Betta, R., Peres, G., Serio, S., & McTiernan, J. 1997, 
A&A, 325, 782 

Ribas, I., Guinan, E. F., Giidel, M., & Audard, M. 2005, ApJ, 622, 
680 

Rocha-Pinto, H. J., Flynn, C., Scalo, J., Hanninen, J., Maciel, 
W. J., & Hensler, G. 2004, A&A, 423, 517 
Rosner, R., Tucker, W. H., & Vaiana, G. S. 1978, ApJ, 220, 643 
Sanz-Forcada, J., Brickhouse, N. S., Dupree, A. K. 2001, ApJ, 554, 
1079 

Scelsi, L., Maggio, A., Peres, G., & Pallavicini, R. 2005, A&A, 432, 
671 

Schmitt, J. H. M. M., Collura, A., Sciortino, S., Vaiana, G. S., 
Hamden, F. R., Jr, Rosner, R. 1990, ApJ, 365, 704 
Schrijver, C. J., Mewe, R., Walter, F. M. 1984, A&A, 138, 258 
Smith, R. K., Brickhouse, N. S., Liedhal, D. A., & Raymond, J. C. 
2001, ApJ, 556, L91 

Soderblom, D. R., Stauffer, J. R., MacGregor, K. B., & Jones, B. F. 
1993, ApJ, 409, 624 

Stauffer, J. R., Caillault, J.-P., Gagne, M., Prosser, C. F., & 
Hartmann, L. W. 1994, ApJS, 91, 625 
Striider, L., et al. 2001, A&A, 365, L18 
Taylor, B. J. 2003, A&A, 398, 731 

Testa, P., Drake J. J., Peres, G., & DeLuca, E. E. 2004, ApJ, 609, 
79 

Tomkin, J., Woolf, V. M., Lambert, D. L., & Lemke, M. 1995, AJ, 
109, 2204 

Turner, M. J. L., et al. 2001, A&A, 365, L27 

van den Oord, G. H. J., Schrijver, C. J., Camphens, M., Mewe, R., 
& Kaastra, J. S. 1997, A&A, 326, 1090 
Vilhu, O. 1984, A&A, 133, 117 
White, N. E., et al. 1994, PASJ, 46, L97 
Withbroe, G. L. 1975, Solar Phys., 45, 301 


Note added in proof. — In a recent publication, Antia & Basu (2005, ApJ, 620, L129) suggested that a significant 
discrepancy between helioseismological studies and modeling of the solar interior, brought about by new measurements 
of a lower solar photospheric abundance of oxygen, may be resolved by an upward revision of the solar photospheric 
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TABLE 1 

Program stars, including a comparison with the Sun 


Star 

Spec. 

type 

Distance^ 

(pc) 

Pvot 

(d) 

log Lx'° 
(erg s-l) 

log Lx'" 
(erg s-l) 

log Lx‘^ 
(erg s-l) 

Age® 

(Gyr) 

47 Gas B 

GO-2 V 

33.56 

^1.0 

30.31 

30.35 

30.39 

0.1 

EK Dra 

GO V 

33.94 

2.75 

29.93 

30.06 

30.08 

0.1 

Tfl UMa 

G1 V 

14.27 

4.7 

29.10 

29.05 

29.06 

0.3 

Ori 

G1 V 

8.66 

5.1 

28.99 

28.95 

28.95 

0.3 

Get 

G5 V 

9.16 

9.2 

28.79 

28.94 

28.95 

0.75 

/3 Com 

GO V 

9.15 

12.4 

28.21 

28.26 

28.26 

1.6 

Sun 

G2 V 

5 X 10-® 

25.4 

27.3 

27.3 

27.3 

4.6 


^stellar distances from Perryman et al. (1997) 
determined from ROSAT in the 0.1-2.4 keV band |Gude^_et_^J1997^[T998a|^ 
^determined from XMM-Newton in the 0.1-2.4 keV band (this work) 
"^dete rmined from XMM-Newton in the 0.1-10 keV band (this work) 
®from[Gude^_et_^y il997ai ll998al^ 


TABLE 2 

Stellar photospheric abundances 


Star 

Ee 

Mg 

Si 

s 

C 

O 

N 

Ref. 

EK Dra 

1.20 

— 

— 

— 

— 

— 

— 

1 

Tfl UMa 

0.83-1.02 

0.65-0.83 

0.78-1.12 

— 

— 

— 

— 

2 


0.93 (0.81-1.07) 

0.74 (0.58-0.95) 

0.89 (0.83-0.95) 

- 

0.85 (0.68-1.07) 

- 

- 

3 


0.87 

— 

— 

— 

— 

— 

— 

1 


1.10 

— 

— 

- 

— 

- 

— 

4 


0.83-0.98 

— 

— 

- 

— 

- 

— 

5 


1.09 (1.00-1.19) 

- 

- 

- 

- 

- 

- 

6 

X^ Ori 

0.89-0.93 

0.91 

0.98 

- 

— 

- 

— 

7 


0.91 

— 

— 

- 

— 

- 

— 

8 


^ 1 

— 

— 

— 

0.63 

— 

— 

9 


1.14 (1.07-1.22) 

- 

- 

- 

- 

- 

- 

6 


1.35 

— 

— 

— 

— 

— 

— 

1 


0.66-1.29 

— 

— 

- 

— 

- 

— 

5 

K,^ Get 

0.89-1.04 

0.85-0.98 

0.95-1.07 

— 

— 

— 

— 

2 


1.29 (1.17-1.41) 
1.0 

1.13 

0.91 (0.76-1.10) 

0.85 (0.78-0.93) 

- 

0.91 (0.71-1.17) 

- 

- 

3 


_ 

_ 

_ 

_ 

_ 

_ 

6 


1.66 

— 

— 

- 

— 

- 

— 

1 


0.98-1.10 

— 

— 

— 

— 

— 

— 

5 

(3 Com 

0.93 

— 

— 

1.38 

0.98 

1.26 (1.05-1.47) 

1.05 

10 


1.00-1.07 

1.17 

1.00 

- 

— 

- 

— 

7 


1.15 (1.07-1.25) 

- 

- 

- 

- 

- 

- 

6 


1.07 

— 

— 

— 

— 

— 

— 

1 


1.00 

— 

— 

- 

— 

- 

— 

8 


1.17 

— 

— 

- 

— 

- 

— 

4 


0.89-1.17 







5 


References. — ■ f2^IOttmann et al.l il99gli : the original values were transformed to solar abun¬ 

dances as giver^^^^^^^^^^^^^^^^g|^^^^^^AG89_exce£^or Fe for which we use the value_^Aen inl^revess^^^Sauv^ll^^Sl 
— GS99; (3) (4)JGra^^i^lJj2TOlJ^; (5) [Ca^re^d^^trobel^^^J (|200lO; f6) lTavlo n2003F !~corrected 

to^^G^99^_fYljE^TOrdasor^r^in il99^) : (8) ; (9) T^ndoir _et_^ !|~|^95 |l corrected to AG89 and GS99; (10) 

IGlegg et al.l 

Note. — Values refer to the solar photospheric composition. If available, error ranges are given in parentheses. 
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TABLE 3 
Observation log 


Star 

Instruments 

Filter 

Start 

Stop 

Exposure [s]^ 

47 Cas B 

RGS 1 

- 

2001/09/10 23:30:24 

2001/09/11 13:39:00 

36610 


RGS 2 

- 

2001/09/10 23:30:24 

2001/09/11 13:39:01 

36610 


MOS 1 

Thick 

2001/09/11 02:04:32 

2001/09/11 13:29:30 

36610 

EK Dra 

RGS 1 

- 

2000/12/30 14:01:58 

2000/12/31 05:17:04 

41960 


RGS 2 

- 

2000/12/30 14:01:58 

2000/12/31 05:17:04 

41960 


MOS 2 

Thick 

2000/12/30 14:10:36 

2000/12/31 04:38:02 

41960 

Tfl UMa 

RGS 1 

- 

2000/11/03 21:44:48 

2000/11/03 12:28:18 

38800 


RGS 2 

- 

2000/11/03 21:44:48 

2000/11/03 12:28:19 

38800 


MOS 1 

Thick 

2000/11/03 21:53:16 

2000/11/03 11:49:04 

38800 

Ori 

RGS 1 

- 

2001/04/07 08:56:49 

2001/04/07 22:31:53 

29326 


RGS 2 

- 

2001/04/07 08:56:49 

2001/04/07 22:31:59 

29326 


MOS 2 

Thick 

2001/04/07 09:03:11 

2001/04/07 17:45:10 

29326 

Get 

RGS 1 

- 

2002/02/09 16:13:01 

2002/02/10 03:21:41 

35920 


RGS 2 

- 

2002/02/09 16:13:01 

2002/02/10 03:21:35 

35920 


MOS 1 

Thick 

2002/02/09 16:19:33 

2002/02/10 03:18:09 

35920 

P Com 

RGS 1 

- 

2003/07/20 02:08:16 

2003/07/20 20:13:28 

61320 


RGS 2 

- 

2003/07/20 02:08:16 

2003/07/20 20:13:34 

61320 


MOS 2 

Thick 

2003/07/20 02:09:06 

2003/07/20 19:12:36 

61320 


^Exposure time used for the analysis, in seconds (excluding flares) 


TABLE 4 

Spectral wavelength 
RANGES USED FOR METHOD 1 


Instrument 

A range (A) 

RGS 

8.30 - 

9.50 

RGS 

12.00 - 

13.95 

RGS 

14.15 - 

15.90 

RGS 

16.20 - 

17.15 

RGS 

17.80 - 

18.30 

RGS 

18.75 - 

19.20 

RGS 

20.80 - 

21.10 

RGS 

21.40 - 

22.40 

RGS 

23.65 - 

24.00 

RGS 

24.50 - 

24.90 

RGS 

28.50 - 

30.10 

RGS 

31.10 - 

32.00 

RGS 

33.40 - 

33.85 

MOS 

1.70'* - 

- 6.90 

MOS 

7.80 - 

9.35 


^For 47 Cas B and EK Dra. 
For the other stars, we used 
lower limits as follows: 4.96 A 
for TT^ UMa and 4.13 A 

for K,^ Cet, and 5.0 A for (3 Com 
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TABLE 5 

Lines used for the EMD reconstruction'^. Method 2, with MEKAL/SPEX 


Line 

logT™^ (K) 

A (A) 

47 Cas B 

EK Dra 

Tfi UMa 

Ori 

Get 

0 Com 

Fe XVII 

6.7 

15.01 

483.2±17.2 

376.9±13.9 

61.5±2.4 

52.0±1.6 

43.5±1.4 

8.5±0.5 

Fe XVII 

6.7 

16.78 

217.5±12.3 

152.7±9.4 

29.3±1.8 

24.1±1.1 

24.8±1.1 

4.3±0.4 

Fe XVIII 

6.8 

14.20 

292.9±15.2 

207.4±11.7 

20.5±1.6 

17.0±1.0 

16.2±1.0 

1.4±0.3 

Fexix 

6.9 

13.52 

222.9±26.6 

126.6±19.6 

14.5±2.9 

8.3±1.6 

6.8±1.5 

0.8±0.6 

Fexx 

7.0 

12.83 

289.7±21.9 

162.5±59.0 

7.4±1.8 

4.3±1.1 

5.8±1.0 

— 

Fexxi 

7.0 

12.29 

210.2±24.7 

149.1±16.3 

8.3±2.1 

7.3±1.3 

8.7±1.4 

- 

Fe XXIII 

7.2 

11.74 

183.3±21.5 

116.5±15.1 

- 

— 

— 

- 

Fe XXIV 

7.3 

10.62 

95.0±40.1 

55.8±13.1 

- 

— 

— 

- 

Fexxv 

7.8 

1.85 

68.0±17.8 

22.6±81.0 

- 

- 

- 

- 

O viii 

6.5 

18.97 

639.8±17.6 

298.5±11.5 

32.2±1.7 

28.2±1.1 

29.5±1.1 

5.6±0.4 

Ovii 

6.3 

21.60 

90.2±11.8 

50.1±8.2 

11.6±1.6 

8.6±1.0 

10.1±0.96 

2.4±0.5 

Ovii 

6.3 

22.10 

48.6±10.2 

39.5±7.7 

8.4±1.5 

6.5±0.8 

6.1±0.9 

3.3±0.5 

Cvi 

6.1 

33.73 

68.2±7.9 

38.1±5.2 

2.6±2.3 

3.6±0.5 

4.1±0.6 

1.3±0.3 

N VII 

6.3 

24.77 

62.2±8.9 

32.3±5.1 

3.8±1.5 

1.8±0.4 

2.5±0.4 

0.4±0.3 

Nex 

6.8 

12.13 

475.0±33.3 

185.7±21.2 

— 

— 

— 

— 

Neix 

6.6 

13.45 

277.8±27.6 

113.0±17.9 

8.4±2.3 

8.9±1.5 

10.0±1.5 

0.9±0.5 

Neix 

6.6 

13.70 

- 

- 

5.9±2.0 

5.2±1.2 

6.4±1.2 

- 

Mg XII 

7.0 

8.42 

238.7±13.2 

93.9±8.1 

3.9±0.7 

2.5±0.4 

6.1±0.5 

0.3±0.2 

Mg XI 

6.8 

9.17 

320.6±10.8 

157.2±7.1 

15.0±0.7 

10.4±0.4 

14.7±0.5 

1.5±0.1 

Si XIV 

7.2 

6.18 

114.2±12.0 

57.3±7.9 

- 

— 

— 

- 

Si XIII 

7.0 

6.65 

194.0±11.0 

105.6±7.5 

5.6±0.6 

4.9±0.4 

5.7±0.4 

0.7±0.2 

S XVI 

7.2 

4.72 

55.6±12.0 

- 

- 

— 

— 

- 

Sxv 

7.2 

5.04 

93.3±11.3 

35.7±6.3 

— 

— 

— 

— 

^For each 

line, the measured luminosity in 10^® 

erg s“^ using 

the MEKAL 

database 

is given. Note 

that the 


entries for the Fe lines contain blends of Fe around the given lines 
^Maximum line formation temperature 


TABLE 6 

Lines used for the EMD reconstruction^. Method 2, with APEC/XSPEC 


Line 

logT^b (K) 

A (A) 

47 Cas B 

EK Dra 

vrl UMa 

Ori 

Get 

0 Com 

Fe VII 

6.7 

15.01 

477.5±19.6 

378.7±16.3 

61.5±2.9 

51.9±1.7 

44.2±1.8 

8.1± 0.7 

Fe VII 

6.7 

16.78 

222.3±18.6 

150.5±15.2 

30.9±2.7 

25.1±1.8 

25.5±1.6 

4.5 ±0.7 

Fe VIII 

6.9 

14.20 

297.4±18.9 

207.1±15.9 

20.4±2.7 

17.2±1.3 

16.9±1.3 

1.4 ±0.6 

Fexix 

6.9 

13.52 

235.2±39.8 

114.7±34.0 

14.4±2.3 

5.8 ±2.8 

6.6±2.2 

0.7±0.6 

Fexx 

7.0 

12.83 

246.8± 30.9 

144.7±30.1 

8.6±1.8 

4.4±1.9 

7.1±1.5 

— 

Fexxi 

7.0 

12.29 

190.0±40.3 

132.2±29.6 

12.3±3.7 

6.8±2.3 

9.7±1.5 

— 

Fexxiii 

7.2 

11.74 

155.4±20.4 

109.6±14.3 

- 

— 

— 

— 

Fe XXIV 

7.3 

10.62 

71.3±59.4 

53.5±33.5 

- 

— 

— 

— 

Fexxv 

7.8 

1.85 

72.7 ±17.9 

40.6±14.1 

- 

- 

- 

- 

O VIII 

6.5 

18.97 

651.1±17.7 

291.4±11.4 

32.9±1.7 

28.4±1.1 

30.1±8.7 

5.6±0.4 

Ovii 

6.3 

21.60 

99.7±20.8 

45.5±13.3 

11.8±1.6 

9.1±1.4 

10.9 ±1.4 

2.8±0.9 

Ovii 

6.3 

22.10 

43.8±10.6 

39.5 ±11.8 

8.4±1.5 

6.4±1.2 

6.5±1.2 

3.1±0.9 

Cvi 

6.1 

33.73 

56.9±7.6 

30.1±5.0 

2.5 ±1.0 

2.8±0.5 

3.5±0.5 

1.3±0.3 

Nvii 

6.3 

24.77 

69.2±9.2 

32.4±5.2 

2.1±0.9 

1.5±0.4 

2.7±0.5 

0.2± 0.2 

Nex 

6.8 

12.13 

500.2±35.3 

167.0±21.7 

— 

— 

— 

— 

Ne IX 

6.6 

13.45 

285.0±42.7 

87.5±36.2 

7.9±4.0 

10.2±2.7 

10.6 ± 2.9 

_C 

Ne IX 

6.6 

13.70 

— 

- 

5.1±2.0 

6.0±1.3 

7.1±1.2 

— 

Mg XII 

7.0 

8.42 

245.9±13.6 

108.9 ±16.1 

3.4 ±0.9 

2.1±0.7 

5.5± 0.8 

0.3±0.2 

Mg XI 

7.0 

9.17 

319.5±11.3 

167.5±13.5 

12.9±1.4 

9.6 ±0.7 

13.2± 0.8 

1.5±0.2 

Si XIV 

7.2 

6.18 

130.6± 24.7 

63.0±11.9 

— 

— 

— 

— 

Si XIII 

7.0 

6.65 

207.3± 75.0 

104.5±14.0 

6.5±1.2 

4.5±0.6 

5.5±0.9 

0.7±0.1 

Sxvi 

7.4 

4.72 

52.9±8.1 

— 

- 

— 

— 

— 

Sxv 

7.2 

5.04 

105.1 ±18.9 

34.1±6.4 






^For each line, the measured luminosity is given in 10^^ erg s ^ using the APEC database. Note that the entries 
for the Fe lines contain blends of Fe around the given lines 
^Maximum line formation temperature 

^For /3 Com none of the Ne lines could be reliably measured by fitting <5 lines using the baseline model from method 
1 and the APEC database within XSPEC 
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TABLE 7 

Mean coronal temperatures'^ and ranges 


Star 

log Lx^ 


SPEX 



APEC 


logT 

log Tmin 

log Tmax 

logT 

log Tmin 

log Tmax 

47 Cas B 

30.39 

7.03 

6.51 

7.60 

7.02 

6.52 

7.65 

EK Dra 

30.08 

6.96 

6.46 

7.45 

6.99 

6.47 

7.58 

ttI UMa 

29.06 

6.65 

6.16 

6.98 

6.69 

6.18 

7.17 

Ori 

28.95 

6.64 

6.14 

7.03 

6.66 

6.22 

6.98 

Get 

28.95 

6.66 

6.18 

7.05 

6.65 

6.13 

7.10 

/3 Com 

28.26 

6.59 

6.15 

6.92 

6.55 

6.07 

6.89 


^log T is the EM-weighted average of log T, Tmin and Tmax are the minimum and maximum 
temperatures that contain 90% of the EM above and below T, respectively 
^Luminosity in erg s”'' determined with XMM-Newton in the 0.1-10 keV band 


TABLE 8 

Reduced foh- the synthetic best-fit 
spectra® 


Star 

2 b 
'^red 

v2 c 
^red 


d.o.f ® 

47 Cas B 

1.22 

1.37 

1.28 

1091 


1.14 

1.35 

— 

1006 

EK Dra 

1.25 

1.34 

1.25 

769 


1.21 

1.29 

— 

671 

ttI UMa 

1.66 

1.51 

— 

415 


1.28 

1.42 

— 

327 

X^ Ori 

1.58 

1.71 

— 

628 


1.28 

1.39 

— 

514 

Get 

1.40 

1.58 

— 

635 


1.35 

1.43 

— 

538 

/3 Com 

1.33 

1.42 

— 

405 


1.24 

1.27 


313 

^The first 

line for 

a given 

star is 

based on 


MEKAL emissivities, the second line on APEC 
emissivities 

^spectrum obtained from method 1 EMD, us¬ 
ing the regions listed in Table 
"^spectrum obtained from method 2 EMD, it¬ 
erated to 'ip^ — 1.0 (based on MEKAL) 
^spectrum obtained from method 2 EMD, it¬ 
erated to '0^=0.5 (based on MEKAL) 

^degrees of freedom 
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TABLE 9 

Abundances with respect to Fe, and absolute Fe abundances'^ 


Abundance (ratio) 

Method 

47 Cas B 

EK Dra 

Tfi UMa 

Ori 

Get 

/9 Com 

C/Fe 

Ml. MEKAL 

0.95t°;« 

0 54+'^' 

QY 

0-291^6 

0 34+0-08 
o’-o^_o Qg 

u.zy _0 lY 


M2. MEKAL 

0.80t°;« 


0 19+0-24 
'0-l^-0.24 

0.34l°+ 

0 43+0-12 
0 '-^o_Q 12 

U.OD _0 29 


Ml. APEC 

0.76t°;« 

U.DU_0 24 

u.^o_Q 13 

i^-iy_o.o5 

n 99 + 0.08 
0 .^^_Q 06 

0 lO^’O-iO 
^•-‘-'-'-0.14 


M2. APEC 



0 14 + 0.09 

0'-^^0.07 

o.^‘*-0.08 

u.oy-o 22 

N/Fe 

Ml. MEKAL 

1 1 Q + C.iiii 

0-84^0.;? 

0 36+0--*^^ 
U.OD_q 12 

0 25+o-o<' 
o.^o_Q Q 6 

0 43+o--<u 

0^-^O_Q QQ 

0 12 + 0 --^^ 
0 '-l^_O .12 


M2. MEKAL 

-0.28 

0.78l°;2« 

0 44+0-28 
^•^^-0.28 

0 26^^’®^ 
Qg 

0 42+0-12 
O'-^^- 0 .i 2 

U.OU _0 29 


Ml. APEC 

-‘-•^'-'- 0.20 

0 71 + 0-26 
^•'-*^- 0.21 

'0-'J^_0.i7 

n 1 c+O-06 

b.iO_Q Q 4 

0 

'^•^‘i-O.oo 

O'-OO _0 05 


M2. APEC 

-^•-^^-0.30 

0 73+0-21 
''^- 0.21 

14 

0 20 ^^’®® 
og 

0 33 + 0.11 

oi-oO-o 11 

0 . 101 °;" 

O/Fe 

Ml. MEKAL 

o.79+";"g 

n 6 n"^o-oy 

U.DU_0.09 

0 27+0-0''' 

'-'•^'-0.06 

0 35+0 05 

u.oo_Q Q 6 

0 36+0-0'' 

U.OD_q Q 6 

0 . 281 °;" 


M2. MEKAL 

0.70l°;ll 

0 51 +0-15 

U.01_o.i5 

n qQ+0-09 

u.oz_0 09 

r\ qq + O.OS 
U.OO_Q Qg 

n 39+0-09 
o>.oy_Q Q 9 

0 41 +0.14 


Ml. APEC 


0 41 + 0-10 
^•^-*^- 0.10 

0 1 7+0.06 
tJ-i ^ -0.09 

0 17+0-03 
^•-*-<-0.05 

n 1 0+0.04 
d.lO-0.04 

U.UD_q og 


M2. APEC 

0.64t°;“, 

n 40+0-07 

U.^U Q QY 

0 25+0-06 
^•^U-0.06 

n oc+0-06 

U.45 q Q 6 

n 9Y+0.09 

oi.^/ -0.09 

n 26'*”^'^^ 

0 '-^D-o.i 2 

Ne/Fe 

Ml. MEKAL 

1 7 Q+c.iy 

1 05+01^ 

0 30+0--*^-*^ 

O'-'^^-O.lO 

0 QO + O-Ob 

n 50+O10 

0 '-OU_o 09 

0 094'O-i^ 

u.uy _0 10 


M2. MEKAL 

1.DO_q 29 

^•'-'-^- 0.21 

0 62"*”®'^^ 
' 0 -o^_ 0 . 2 l 

Q 70 + O.I 6 

O'- ''J_0.16 

u.yo_Q 20 

0 41+0.35 

'0-^i-0.35 


Ml APEC 

1 

-‘-•'^^-0.18 

U.yO_Q 25 


rv 09 + 0.09 

u.oz_0 07 

0 44+0-10 

'O-^^-O.lO 

'o-oo_Q 10 


M2. APEC 

1 yc+O.is 

0 79+0-17 
^•'^-0.17 


0.59l«;i= 

o'-o)y o 19 

- 

Mg/Fe 

Ml. MEKAL 

1 yi 4 -U.lY 

1 oi+U.i/ 
-L-oI-0.17 

1 OO + '-^-^y 
-^•^^-0.27 

0.981“;" 

1 49 +o-^o 
-'-•'*^-0.23 

1 00 '*’O-^^ 
-^•oo'_o.5l 


M2. MEKAL 

r) Q -1 +0.38 
^•'^-L_o.38 

1 54+0-28 
l-O^_0.28 

1 24"*”®'^^ 

1-^4_q 31 

1 1 9+0.18 
i-i^-0.18 

1 94+0-80 
i-y4_0 30 

i. 2 ol 2 ;?l 


Ml APEC 

J--OU_o 16 

1 oi +0.26 

^•^-^-0.29 

1 03+0-53 

l.UO_Q 54 

n 5^1 + 0-21 
^•°-*^- 0.20 

1 17+0-27 
' -0.27 

0.901^2 


M2. APEC 

1 97+0-23 
-0.23 

i-4ho:26 

0.73t°;2i 

0 73+0.17 

O'- <0-0 17 

l-l+n.ls 

n 81 +0-18 
b-oi-o 18 

Si/Fe 

Ml. MEKAL 

o.89+^;li 

o.89t":;^ 

0.69+^;i? 

0.74l“:" 

0.871“;" 

q+U.b4 

-'-•J-O- 0.68 


M2. MEKAL 

i-07tg;h 

0.95t°:i| 

0 . 68 l°;“ 

0.82l°:l® 

o.95l°;l? 

n 95 + 0-88 

u.yo_o 38 


Ml APEC 

n 79+010 

0 >- ' 1 ^_ 0.10 

0.85l°;« 

0.60l°i^ 

o.eel®;!® 

o.73lg;h 

0 77+0.63 


M2. APEC 

1 .UD o .20 

o-7»-o:\l 

0.50l°;“ 

0 57+0-11 
0 '- 0 /o.ii 

0.63l°;“i 

0.69l°'J7 

S/Fe 

Ml. MEKAL 

0 77+^--*-^ 

‘ ‘ -0.14 

0 59+^-^^ 

u.oy_Q 15 

- 

- 

- 

- 


M2. MEKAL 

±.uo_o 32 

0 59+0-18 

u.oy_Q ig 

- 

- 

- 

- 


Ml. APEC 

'^•'^'-'- 0.11 

0.54l°;« 

- 

- 

- 

- 


M2. APEC 

-^•*^^-0.26 

0.50t°:i® 

- 

- 

- 

- 

Ar/Fe 

Ml. MEKAL 

1 c;o+0.43 

i.OO_Q 42 

0 44+0-44 
0.44:_o 44 

- 

- 

- 

- 


Ml. APEC 

-L-^y-0.33 

n 55+0-47 
^•O0 o.45 

- 

- 

- 

- 

Fe 

Ml. MEKAL 

n 51 +0-03 

'J-oi-0.04 

n 63+0 00 

O.Dc5_o.o5 

o.73+t;;M 

o.63l“;“® 

Q 71 +U.US 

OJ-/1_Q.06 

0.431°;" 


M2. MEKAL 

OI-OU_0.05 

r\ yo + O.OS 
^•'^_0.08 

0 . 8 ll°;h 

0.87l°:h 

i.isig'it 

1.27ll;g 


Ml. APEC 

n 69+0-00 
U.Dy_0.05 

o-74t;^:J^ 

1 07+°-^® 

1-Of _0.32 

o-98lg:h 

o.93l°;l? 

1.7612;!? 


M2. APEC 

n 55+0-05 
^•oO-0.05 

0.96l°+ 

1 26'*”^'^^ 
1 -^ 0 _Q 31 

0.83l®;i| 

1 qq+0.48 

1 . 00 _Q_ 4 g 

— 


^All abundance ratios and Fe abundances are with respect to the solar photos£heric_abundances ei^n by 
|Ander^_&_Greves^ (198^ except for Fe, for which the photospheric value given by |o revesse_&_^auv^ (199^ has 
been adopted 
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TABLE 10 

Light curve modeling 



47 Cas B 

EK Dra 

ttI UMa 

Ori 

Cet 

f3 Com 

Slope for T < 

2.10-2.63 

2.22-3.09 

1.51-2.46 

1.88-2.86 

1.85-2.80 

2.12-3.67 

Slope for 

(-1.34)-(-1.48) 

(-1.55)-(-1.62) 

(-2.45)-(-3.63) 

(-3.26)-(-3.68) 

(-2.25)-(-2.60) 

(-2.95)-(-3.53) 

4> 

-0.3 

-0.3 

-0.3 

-0.3 

-0.3 

-0.3 

C 

0.95-0.76 

0.90-0.65 

1.32-1.23 

1.06-0.70 

1.08-0.71 

0.94-0.54 

a{l3 = 0) 

2.21-2.28 

2.25-2.31 

2.43-2.78 

2.60-2.80 

2.39-2.54 

2.54-2.76 

a(/3 = 0.25) 

2.16-2.20 

2.19-2.23 

2.33-2.59 

2.45-2.60 

2.30-2.41 

2.40-2.57 

\ogE2/Ei{l3 = 0) 

(-4.2)-(-5.3) 

(-2.8)-(-3.4) 

(-1.8)-(-2.2) 

- 

- 

- 

log£;2/Ei(/3 = 0.25) 

(-6.2)-(-6.8) 

(-3.5)-(-4.1) 

(-1.8)-(-2.5) 

— 

- 

- 

log E 2 

30.20 

30.00 

29.0 

— 

— 

— 

mod. depth^ 

0.028 

0.056 

0.122 





^Determined in the logarithmic energy range between logT = 6.2 and (log Tm-O.l) 

^Determined in the logarithmic energy range between (logT^+O.!) and an upper limit that depends on the star (logT = 7.5 for 
47 Cas and EK Dra, logT — 7.3 for tt^ UMa, and logT — 7.1 for the other stars) 

^Modulation depth of the observed light curves 
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Photon energy (keV) 


Fig. 2.— MOS spectra of the six solar analogs, normalized to a distance of 33.56 pc, the distance of 47 Cas B. The overall luminosity of 
the corona decreases from top to bottom as indicated by the labels. 
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Fig. 3.— Upper panel: data and fitted spectrum for 47 Cas B using method 1 with APEC; lower panel: same for method 2. The 
synthesized spectra from the best-fit parameters are shown in black for the ranges that were used in method 1. Error bars are displayed 
only for the MOS data in order to avoid confusion. The insets show the important spectral portion containing the Fe L shell lines. The 
lower plots in each figure show the contributions of each bin to the value. Note that for method 2, we show the identical layout for 
illustration and comparison purposes although the iteration procedure did not make use of a binned spectrum. 













































Emissivity [10“ erg cm^ /sec] 


26 


Telleschi et al. 


Emissivities of Fe and O 



6.0 6.5 7.0 7.5 8.0 

log(T emperature) 


Fig. 4.— Emissivities f or Fe and O line blends from the MEKAL database, plotted on a logarithmic grid, assuming solar abundances 
lAnders fc Grevesselll98tl) . 
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Fig. 5. — Reconstructed emission measure distributions. Left: Prom method 1 using SPEX, based on Chebychev polynomials of 
order 6 (solid line) and order 8 (dashed). Middle: From method 2, based on MEKAL emissivities. Right: From method 2, based on 
APEC emissivities. In the middle and right plots, the red histograms illustrate the ± Icr range of solutions from the average of 20 EMDs 
reconstructed from the original and from the perturbed line lists. The black histograms illustrate the best-fit EMDs, derived from the 
unperturbed line-flux list. 
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Fig. 6. — Fractional deviation of the predicted Fe line fluxes and the O vill/O VII flux ratio from the observed values for method 2. Filled 
circles: SPEX. Open circles: APEC. 
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Fig. 7.— Emission measure distributions (using SPEX with method 2). Left: The iteration was stopped when = 1.0 was reached. 
Right: The iteration was stopped when 'ip'^ = 0.5 was reached. Note the additional oscillations in the EMD. 
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Fig. 8.— Abundan ces relative to Fe as a function of FIP, normalized to solar photospheric ratios lAnders &: Grevessd 119891 

IGrevesse &: Sauvallll993l . Open circles: method 1; filled circles: method 2. 
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Fig. 9. — Light curves of our targets. The four light curves in each panel show, from top to bottom, the total count rate in the 0.2 — 10 keV 
range (black), in the soft band (0.2 — 1 keV, green), in the hard band (> 1 keV, red, where the upper energy limit is reported in the caption 
of TableUI, and the ratio of hard/soft (blue). For illustration purposes, the hardness ratio has been multiplied by 5, 3, 0.5, 0.5, 0.5, and 0.5 
for 47 Cas B, EK Dra, tt^ UMa, Get, and {3 Com, respectively. The bin size is, for the stars as listed above, 300, 300, 300, 450, 

450, and 600 s, respectively. Only data from detectors that were operated in imaging mode were used, i.e., data from the PN camera were 
not used for tt^ UMa, Cet. For the latter two stars, only one MOS camera was available in imaging mode (see Table 

The time ranges of the largest flares excluded from the spectral analysis are also shown by dash-dotted vertical lines. 
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Fig. 13. — Abundances of Fe, Ne, and abundance ratios of Ne/Fe, O/Ne, O/Fe, and Mg/Fe are plotted as a function of the mean coronal 
temper ature. Value s from method 2 (MEKAL) have been used. The dotted contours delimit the regions derived from a larger stellar 
sample l(;udeJI2004ll . Note that the displayed range of the abundances or abundance ratios is 1 dex in each of the six plots. 
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71' UMa 




Fig. 14. — Simulated light curves obtained assuming a flare distribution that is based on the a values found fro m our EMDs. The 
maximum flare energy and the difference between largest and smallest flares assumed here are given in Table irni The black shapes 
represent the largest flares actually used to synthesize the light curves from the simulations (and this flare may therefore be smaller than 
the largest flare actually observed in the light curves in Fig. [9] which are not considered for this comparison). 


















